Ported Tutte Functions of Extensors and Oriented 

Matroids 

Seth Chaiken * 

Computer Science Department 
University at Albany, State Univ. of New York, USA 
sdcOcs . albany . edu 

August, 2 2006 

Mathematics Subject Classifications: 52C40, 15A75, 05C50 



*Some of this research was performed during my 200f academic year Sabbatical at the University at 
Albany. 



1 



Abstract 



The Tutte equations are ported (or set-pointed) when the equations F(J\f) = 
g e F(N/e) + r e F(N\ e) are omitted for elements e in a distinguished set called 
ports. The solutions F, called ported Tutte functions, can distinguish different 
orientations of the same matroid. A ported extensor with ground set is a (fully) 
decomposable element in the exterior algebra (of antisymmetric tensors) over a 
vector space with a given basis, called the ground set, containing a distinguished 
subset called ports. A ported extensor is one way to present a linearly representable 
ported matroid or oriented matroid. There are extensor operations corresponding 
to oriented matroid dualization, and to deletions and contractions. 

We define a ported extensor function by means of dualization, port element 
renaming, exterior multiplication, and then contraction of all non-port elements. 
The main result is that this function satisfies a sign-corrected variant of the Tutte 
equations in which deletion and contraction are extensor operations, and addition 
and the anticommutative multiplication belong to an exterior algebra rather than 
to a commutative ring. 

For extensors representing unimodular, i.e., regular matroids with an empty port 
set, our function reduces to the basis generating function; and then, for graphs, to 
the Laplacian (or Kirchhoff) determinant. On graphs with port edges, the function 
value, as an extensor, signifies the space of solutions to Kirchhoff's and Ohm's 
electricity equations after projection to the voltage and current variables associated 
to the ports. In particular, the Laplacian matrix with the identity matrix appended 
presents one example of our extensor function's value. Combinatorial interpretation 
of various determinants (the Plucker coordinates) generalize the matrix tree theorem 
and forest enumeration expressions for electrical resistance. 

We also demonstrate how the corank-nullity polynomial, basis expansions with 
activities, and a geometric lattice expansion generalize to ported Tutte functions of 
oriented matroids. The ported Tutte functions are parametrized, which raises the 
problem of how to generalize known characterizations of parametrized non-ported 
Tutte functions. 
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1 Introduction 

A ported unoriented or oriented matroid M = J\f(P, E) has its ground set S(Af) = P U E 
given with a distinguished subset P of elements which we call ports; P D E — 0. The 
following definition combines the idea of Tutte invariants of set-pointed matroids studied 
by Las Vergnas [35-37] with the idea of the parametrized Tutte equations and functions 
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studied by Zaslavsky [67] (who calls our Tutte functions "strong"), Bollobas and Riordan 
[5] and Ellis-Monaghan and Traldi [23]. Let two parameters g e and r e be given for each 
e e E. 

Definition 1.1. A function F is a (ported and parametrized) Tutte function if 

the domain of F is a minor closed class of ported unoriented or oriented matroids and 
F satisfies the following (ported and parametrized) Tutte Equations for each M in 
the class: 

When e G E is a non-separating element, i.e., e is neither a port nor loop nor a coloop 
(i.e., isthmus): 

F{M) = g e F{M/e) + r e F{M\e). (1) 
When Af = A/i © Af 2 with ground sets ^(A/i) H S(JV 2 ) = 0: 

F(Ni © M 2 ) = F(A/"i)F(A/" 2 ). (2) 

In the following, a Tutte function shall be assumed to be both ported and parametrized 
as above unless otherwise indicated. The letters g e and r e which we use for parameters, and 
the terminology "set of ports" for the distinguished ground set elements, are suggested 
by the electrical network theory application which motivates our research. Briefly, g e 
(for conductance or admittance) and r e (for resistance or impedance 1 ) appear in the 
homogeneous expression r e : g e of a resistance value (i.e., the resistance value is either 
r e / g e ohms if g e ^ 0, or infinity if g e = and r e ^ 0.) A port designates one interface (a 
pair of terminal vertices and two variables, one for current and one for voltage) between 
the physical electrical network and an external environment. (Think of an ordinary 2-pin 
electrical plug or receptacle.) Other physical applications of parametrized Tutte functions 
of graphs are surveyed in [48]. Details and references for the electrical application are 
presented but the logical background, results and proofs (^210) are independent of these 
details. 

Our main topic is a generalization of one construction of a linear subspace from the 
graph (with port edges and parameters) that occurs in the application. The constructed 
subspace consists of the space of network solutions projected to the port variables. We 
generalize by replacing the graph's coboundary (cocycle) space by a finite dimensional 
linear subspace. The underlying space has a particular ground set basis, whose elements 
generalize graph edges, and this ground set has the distinguished subset of port elements. 

Our formulation with ports subsumes (see the classical equation L(g, 1)0 = J 

on the node voltages {0.;} and external currents flowing into nodes { Ji}. Here, L(g, 1) 
denotes the edge- weighted Laplacian matrix; the edge weights g e are the conductances (or 
admittances) of edges and we take each r e = 1. Assuming the graph is connected, each 
principal cofactor of L(g, 1) equals (according to the famous Matrix Tree Theorem) the 
weighted spanning tree polynomial 

t( 9 )= e n*=i>- 

TCE eST T 

T a spanning tree 

1 Electrical engineers customarily use the words admittance and impedance for complex values. 
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Choe, Oxley, Sokal and Wagner [66] prove that this and similar basis generating poly- 
nomials for some but not all other matroids have the half-plane property: For complex 
values g e , if 8lg e > for all e 6 E then T(g) ^ 0. The half-plane property charac- 
terizes those electrical networks that cannot generate energy. These authors show that 
the monomials of every non-zero homogeneous multiamne polynomial with the half-plane 
property enumerate the bases of a matroid. After noting that the half-plane property is 
a strengthening of the property that all the coefficients have the same phase, they raise 
the open question of whether, for every matroid, coefficient values exist so the basis enu- 
merating polynomial with these coefficients has the half-plane property. It is noted that 
such coefficients can always be assumed to be positive. 

We treat polynomials (including T(g) when r e = 1) whose ratios, for graphs, signify 
externally observable numerical characteristics of an electrical network that pertain to sev- 
eral port current or voltage values together. These quantities include the "transpedances" 
used by Brooks, Smith, Stone and Tutte [7] in their work on finding dissections of squares 
into squares by means of solving electrical network equations. Indeed, these authors' 
presentation of the Matrix Tree Theorem is in the context of combinatorial proofs of 
solution properties. Our treatment describes the coefficient signs in terms of oriented 
matroids. Some of the polynomials have the half-plane property and some clearly do 
not because they have terms with different signs. We therefore hope that our work may 
contribute insight into the questions posed in [66] , as well as more recent questions raised 
by Wagner [58] pertaining to inequalities among differences between polynomials that 
enumerate certain trees or forests. This work combines Rayleigh's inequality (motivated 
by the physics of electrical networks with positive edge conductances, see also [17]) and 
analysis of correlations where ratios of edge parameters signify probabilities. 

It may help some readers to know that the generalization of the coboundary space and 
the parameters r e and g e , which together determine our constructed solution subspace, are 
specified separately — This separation corresponds to the distinction between two kinds of 
constraints (exact or geometric versus approximate) that is explained and modeled with 
matroid theory by Murota [43]. 

Our work distinguishes the polynomial determinants where all terms (in the param- 
eters) have the same sign from those where differing signs occur. So, a non-zero deter- 
minant with terms of differing signs might vanish for critical combinations of parameter 
values; whereas a determinant with terms all the same sign will never vanish. In Murota's 
model, this distinction would depend on the exact or geometric constraints; we express 
the distinction in terms of its oriented matroid properties. 

Some of our previous work [13] applied oriented matroids to distinguish the case of a 
vanishing determinant for critical parameter value combinations. The approach did not 
apply Tutte function theory. Instead, we investigated the property of a pair of oriented 
matroids with a common ground set have a common non-zero covector. The electrical 
network applications treated were more general than in the present paper. We began 
with a dual oriented matroid pair (graphic and cographic), deleted or contracted certain 
elements in each, and then evaluated the above common covector property in the resulting 
oriented matroid pair with a common ground set. See [62] for Tutte theory developed 
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for paired matroids. We are currently investigating how to generalize paired matroids by 
adding port elements and relate the theories of the present paper with [13] and [62]. 

1.1 Exterior Algebra 

Exterior algebra is used to represent and operate on linear subspaces. A (fully or com- 
pletely) decomposable element in an exterior algebra is either a field element or the exterior 
product of vectors. For the sake of brevity, such decomposable will be called extensors. 
An extensor is ported when an underlying space is given with a ground set basis and a 
distinguished port subset. We remind the reader that each non-zero extensor corresponds 
to the unique subspace whose bases are the sets of vectors whose exterior product equals 
the given extensor, up to a non-zero field element multiple (see Theorem 12.11 ) 

For the electrical network application, a graph's coboundary space is represented by 
extensor N. In our theory, N would represent an arbitrary linear subspace. We will 
define some functions of N or of certain matroids with values within an exterior algebra. 
Each function value represents the solution subspace for a generalization of the electrical 
network problem. The domain of the first such function (see Definition IT2"|) consists of 
ported extensors. The second function applies to the class of ported unimodular (i.e., 
regular) matroids. The second function is defined by specializing the first function to 
the extensors that represent unimodular matroids. See Definition 14.31 When the latter 
function is specialized to graphic oriented matroids with no port elements, it reduces to 
the Laplacian determinant which, according to the famous Matrix Tree Theorem, equals 
the spanning tree polynomial discussed above. 

The additive Tutte equation pertains to a graph or matroid and its two minors obtained 
by deletion and contraction. We will define algebraic operations on N. Their values are 
denoted by N/e and N\e for each e G E. We also adapt the Hodge star operation to define 
an operation corresponding to oriented matroid dualization. Our main results (Theorems 
13.31 and 14 .5|) are that these two functions obey equations similar to the Tutte equations 
above, taken in an exterior algebra rather than in a commutative ring. The dualization 
operation is used in defining these functions and in proving the main results. We hope the 
reader will bear with us in using the terms Tutte equations and functions in this context 
before the precise equations can be presented. If not, one can skip to the definitions and 
theorems in and ^21 without loss of logical continuity. Sign factors are required in our 
Tutte equation variants to accommodate the anticommutative multiplication. 

It is well-known and easy to verify that the spanning tree count is a Tutte invariant 
of graphs. Our results further elucidate the relationships between Tutte functions, the 
Matrix Tree Theorem theorem and enumeration methods for resistive electrical network 
solutions pioneered by Kirchhoff [33] and Maxwell [42]. These methods were introduced 
into combinatorial theory by Brooks, Smith, Stone and Tutte [7,53] who attributed them 
to Kirchhoff, see ^6.11 They continue to be applied within some electrical engineering 
computer aided design tools [24]. Our two points of departure from [7] are to replace 
analysis in terms of graph vertices and incidences by analysis of functions on the graph 
edges, and then to express the relevant equations f ^6.2|) in exterior algebra. 
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The extensor functions that we study are generalizations of the determinant of the 
(reduced) edge-weighted Laplacian matrix of a graph. The homogeneous form of this 
determinant, according to the Matrix Tree Theorem, equals 

r s &et(NDN t )=r s Y J 9 ^ L N[X? = rsYs 9 f = Y.^ r B> ^ = diag(^, e G £) (3) 

XCS TX B TB B Te 

where N is the reduced signed incidence matrix of the graph and the sum is over all 
spanning spanning trees B. We remind the reader that every extensor can be represented 
by a sequence of determinants which is called its Pliicker coordinates. Of course, the 
maximal forests are the bases of the graphic matroid. 

In summary, our theorems about extensor functions satisfying sign-corrected Tutte 
equations generalize the graphic matroid case of the easy-to-prove fact that the basis 
generating function is a (non-ported) Tutte function defined for all matroids. 

1.2 New Tutte-Like Invariants 

However, somewhat deeper theory involving the Laplacian, matroid invariants and exte- 
rior algebra is involved. First, consider Tutte invariants of matroids, such as the number 
of bases. The universal Tutte invariant, the well-known Tutte polynomial, is defined for 
all matroids; hence no Tutte invariant depends on the orientation of an oriented matroid. 
Let us extend the definition as follows: A ported Tutte invariant of ported oriented 
matroids is a Tutte function / where all the parameters r e = g e = 1 that is invariant 
under oriented matroid isomorphisms / that preserve each port element. Specifically, 
f(p) = / _1 (p) = P f° r every port p in the domain or range of /. Las Vergnas [35-37] 
developed the analogous ideas for ported matroids (without orientation), and called them 
"set-pointed." This work on invariants of matroid morphisms, i.e., strong maps, applies 
the universal ported invariant called the big Tutte polynomial. We had applied this 
idea to study the behavior of such polynomials under the restricted matroid union oper- 
ation [12] where the composed matroids can only have port elements in common, which 
is one generalization of matroid series connection. 

It is easy to develop the corresponding universal ported oriented matroid invariant. 
Since the latter has one variable for each connected oriented matroid whose ground set 
contains port elements exclusively, we conclude that some ported Tutte invariants of 
oriented matroids do distinguish different orientations of the same matroid. We mention 
that J\f is unoriented or oriented in the definition of a ported Tutte function because the 
big Tutte polynomial never distinguishes different orientations of the same matroid but 
the oriented variant of it, which is defined for ported oriented matroids, does distinguish 
some orientations of the same matroid. 

Here is a simple example: Let the set of ports be P = {pi,P2}- Let A/"f (P, 0) be the 
oriented matroid with ground set P and oriented circuit collection {++, — }. A/"f is one 
of the two orientations of the rank 1 uniform matroid with ground set P. Let Af^~(P, 0) 
be the other orientation — Its oriented circuit collection is {H — , — h}. Since and M{ 
are decomposable under ported Tutte decomposition, any ported Tutte invariant F of 
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oriented matroids for which F(Af^~) 7^ F(J\ff) distinguishes different orientations of the 
same matroid. More interesting examples are given in tyj.'Sl 

The basic theory of non-ported matroid Tutte invariants proves that the Tutte poly- 
nomial has various combinatorial interpretations, i.e., different expansions over subsets or 
other structures. The generalizations of these for ported Tutte functions and invariants 
are discussed in — They all include variables identified with (connected) matroids on 
port elements as in Las Vergnas' big Tutte polynomial, except that these matroids are 
oriented. One such expansion generalizes the corank-nullity polynomial. One of our re- 
sults is that our extensor-valued Tutte function of ported unimodular oriented matroids 
can be expressed by replacing each oriented matroid monomial by this function's value, 
an extensor, on the oriented matroid which the monomial identifies. Since the expression 
also requires u — v — where u, v are the corank-nullity polynomial's variables, we see 
that our extensor-valued Tutte function generalizes the basis enumerator. 

Our extensor-valued Tutte function provides one example of a ported Tutte invariant 
of graphic oriented matroids (and more generally, unimodular oriented matroids). Given a 
linear representation N of a graphic or other unimodular oriented matroid, the invariant's 
value is defined by our extensor function applied to the extensor presenting N, when we 
take all the parameters r e — g e — 1. We will see that this generalizes the fact that the 
reduced Laplacian determinant equals the number of bases in the corresponding graphic 
matroid, which is the number of spanning trees when the graph is connected. 

1.3 Structure of Solutions 

The second involvement between Laplacians, ported Tutte functions and exterior alge- 
bra beyond the matrix tree theorem grows out from the discrete (i.e., matrix) variant of 
Laplace's equation, the problems in classical physics that it models, and the structure and 
solution methods for those problems. Among many other analogs, the parametrized, dis- 
crete equation of Laplace combines Kirchhoff's current law (of flow conservation), Kirch- 
hoff 's voltage law (electromotive force in direct current electrical systems is determined 
by a potential function), and Ohm's law (current flow is proportional to potential differ- 
ence). Given a resistive electrical network, the voltages (i.e., potentials) at all vertices 
can be determined, by solving Laplace's equation, from given voltages at certain vertices 
and/or the currents into other vertices flowing from the environment. Our work might 
contribute a few new insights of matroid theory into this situation. One of our starting 
points for developing this application is well-known in both matroid theory and electrical 
network theory: Certain insights are obtained when one can replace analysis in terms of 
graph vertices by analyses involving graph edges and matroids on them. Hence, we model 
the interactions with the environment by currents through, and voltage drops across port 
edges, instead of currents into vertices and potential values at vertices. 

Consider all electrical networks with a fixed set of port edges P. Consider two networks 
to be the same members of this class if they have the same edges and the same cycle spaces. 
Such equivalent networks will have the same electrical behavior. For us, this class is the 
same as the class of graphic oriented matroids with a fixed set of port edges for which 
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two parameters, r e and g e are given for each non-port edge. Each network determines the 
ported extensor that represents the graph's graphic oriented matroid. When our extensor 
function is applied to such an extensor, the function value is an extensor corresponding to 
the linear subspace of solutions projected to the port variables. For the electrical network 
application, the significance of our main theorem is that the function of networks with 
port edges P that gives each network's projected solution space is a ported Tutte function. 
Routine network analysis to solve for port behavior can use the Laplacian determinant 
and cofactors (or equivalent elimination methods) to solve for the constraints between 
port variables. But our results show that not only is the Laplacian determinant a Tutte 
function, but when an extensor is used to express the solution space, the entire solution 
is a Tutte function. As a corollary, the network solution when expressed in exterior 
algebra can be written as linear combination of solutions of networks with no resistor 
edges at all — only port edges. The coefficients are homogeneous multilinear products of 
the parameters. 

The extensor value is can be calculated in three ways. One is by recursive application 
of the Tutte equations, i.e., Tutte decomposition. The second is by substitution into 
a polynomial that extends the parametrized Tutte polynomial to ported, parametrized 
oriented matroids. The third is by Gaussian elimination. That calculation generalizes the 
evaluation of the Laplacian determinant. For graphs it is equivalent to finding a suitable 
representation of the solution space of an electrical network projected onto the variables 
associated to the port edges. Of course, our extensor-valued invariant is calculated after 
setting all r e — g e — 1. 

1.4 Computational Complexity 

An additional motivation for our work comes from theoretical computer science. The 
number of bases in unimodular matroids is virtually the only non-trivial Tutte invariant 
whose computation is tractable. For such matroids, the computation input is a totally 
unimodular matrix N and the output value is det NN*. The more general problem, 
to compute the number of bases a matroid represented by an arbitrary matrix, is j^P- 
complete (see $ 17.21 for details.) Our generalization of this detiViV* is an extensor; a 
succinct matrix representation of it can be computed by Gaussian elimination (like a 
determinant). Therefore, when the notion of Tutte invariant is generalized with ports and 
with exterior algebraic values (and the Tutte equations are modified with sign corrections), 
new computationally tractable invariants are obtained. 

1.5 Additional Context 

See [43, 46] for elements of electrical network theory from the point of view of matroids. 
Electrical networks with non-linear but monotone resistance functions were studied using 
graphs by [28,29] in a way that led us to apply oriented matroids to this topic [13]; these 
publications may help orient the reader our point of view. Our current work and [12] 
extend to sets of more than one distinguished element results and ideas about series 
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and parallel connection of both networks and matroids. In particular, the four-variable 
Tutte polynomial for a pointed matroid, which was defined and applied to these topics 
by Brylawski [8], is a special case of our Tutte function. Our previous work covered 
solutions [11] by energy minimization [39] of electrical networks (and generalizations) with 
multiple ports, applications of enumeration [15] combined with oriented matroid ideas [14] 
to network analysis, and a (non-oriented) matroid abstraction [10] of the solutions. 

We proceed to the details about exterior algebra and Pliicker coordinates pertaining 
to realizable non-oriented and oriented matroids. They include algebraic operations and 
identities which correspond to some elementary matroid relationships. Our main con- 
struction, an extensor valued function of ported extensors, and the ported Tutte equation 
variant that it satisfies, is presented in A variant of the corank-nullity polynomial in 
§Uis used to express our function restricted to unimodular extensors. The variant differs 
from Las Vergnas' big Tutte polynomial so (1) it applies to oriented matroids instead 
of unoriented matroids, and (2) it includes parameters as in Definition 11.11 Extensions 
to ported matroids of known results about expressing Tutte functions as set, basis and 
flat expansions, and the related open questions follow in £0 Further discussion of the 
electrical context is given in §U1 and brief remarks on peripheral topics appear in ^7.31 

2 Preliminaries 

Throughout, K denotes a field, either the reals, rationals, or their extensions generated 
by the parameters g e , r e . 

2.1 Exterior Algebra 

We refer the reader to basic texts such as [31, §7.1-7.2, on associative and exterior alge- 
bras over fields] for complete development and proofs. The following is a synopsis with 
the emphasis on the facts we will need. It also explains certain notational conventions 
which help to mimic oriented matroid theory in exterior algebra. We use a combinatorial 
approach to adapt the operations of Hodge star (for duality) and tensor contraction (for 
matroid contraction). 

An associative algebra A over K is a ring that is also a vector space over K, for which 
addition and are the same in both the ring and the vector space, and for which the 
ring and scalar multiplications are compatible: a(xy) = (ax)y = x(ay) where a G K and 
x, y G A. 

Let V be the vector space KS where finite set S is a basis. Thus V consists of the all 
X] e es aee ' ° e ^ where a e e = if and only if a e = for every e G S. The associative 
algebra over K generated by V consists of all finite if-linear combinations of 1 (the ring 
identity) and formal finite (non- commuting) products of elements of S. 

The exterior algebra £(V) over V is the quotient of the associative algebra over K 
generated by V modulo the algebra ideal / generated by products v 2 , v G V . The image 
of each non- zero v G V under the map v — > v + / G S(V) is denoted by v. These v 
will also be called vectors. Thus, for Vi,v% G V, (vi + v 2 )(vi + v 2 ) = 0, VjVj = and 
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so V!V 2 = — v 2 V! in £(V). This and the associativity law imply that each product of a 
sequence of vectors vanishes if the sequence has repeated elements. Indeed the product 
vanishes if and only if there is a linear dependency among the vectors being multiplied. 
Note that a non-zero product of two or more vectors is not a vector. 

A particular basis of £(V) is constructed from the basis S = {si, . . . , s n } of V. This 
basis consists of 1 together with the 2 n — 1 products of vectors in distinct non-empty 
subsets of {si, . . . , s n }, each product written in a particular order. A formula for exterior 
product expressed in terms of this basis is used in [31] to prove that the product is 
associative. The formula expresses the following fact for products of basis vectors which 
is true for all products: Given any sequence of vectors Vi ■ ■ ■ and permutation cr £ & k 
with sign e(er), the exterior product satisfies the alternating law 

vi ■ • • v fc = e(cr)v ai ■ ■ ■ v CTfc . 

As a result, £(V) is an associative algebra that has dimension 2 n when viewed as a vector 
space over K. It is customary to use increasing order of subscripts, so each X £ £(V) 
can be expressed by 

X = xqI + ^2 x A s h ■ ■ ■ s i k 

%^ACS: 
A={s il ,...,Si k }, h<—<ik 

with 2" unique coefficients xa, A C S. 

We follow a different convention which mimics the one used with the chirotope crypto- 
morphism for oriented matroids given in [4] . For us, A will denote an arbitrary sequence 
of elements of S. The value in K symbolized by coefficient xa will depend on the order 
as well as the elements of A, but these values will satisfy 

XA a = e(a)x A 

where A a = (oi • • • a k ) a = a a i • • • o-ak is A permuted by a. In general, a function like 
A — ► xa is called alternating if it has this property. Our convention allows A to have 
repeated elements but the alternating property requires x A = for such A. 

We follow a related convention for subset expansions and formulas within them. When 
necessary, a set symbol like A in A C S will denote distinct elements written in an 
arbitrary sequence. But the expansion or formula will be written only if its value is 
independent of the sequence chosen for each symbol. Furthermore, when A is a sequence of 
distinct basis vectors, the corresponding product of their images in £(V) will be denoted by 
A. The empty sequence corresponds to 1 £ £(V). No sequence of linearly independent 
vectors corresponds to £ £{V). The above 2 n term basis expansion is thus written 

X= ^ X A&1 ■ ■ ■ SL\A\ = ^2 XaA - 

ACS ACS 
A={ai,...,o| A |} 

This expansion follows our convention because xa ct A ct = 6 2 (<j)xaA- = x^A. Note that 
X = if and only if every x A — 0. 

The concatenation of sequences A, B, C, . . . is denoted by ABC . . . 
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2.2 Extensors with Ground Set 



The exterior algebra £(V) is a powerful tool to explore, in a coordinate free way, relation- 
ships between linear subspaces of V [2,65]. These relationships are, in other words, the 
theorems about the projective geometry whose points are the rank 1 (zero-dimensional) 
subspaces. The geometric flats (the empty set, points, lines, planes, etc., and the whole 
space) correspond to these linear subspaces. One way to present a i^-realizable matroid 
is to map each ground set element to either a point or to the empty flat in this projective 
geometry. The matroid structure is then expressed in terms of incidence of the images of 
the ground set elements with the geometric flats. The formulation is coordinate-free be- 
cause these relationships do not change under a change of the basis S for V. We mention 
this formulation to contrast it with our application of exterior algebra. 

In our application, each i^-realizable matroid with ground set S will be presented by 
a separate (fully decomposable) non-zero element N £ £(V), where V = KS. Each such 
element N will determine a linear subspace L = L(N) of KS. Consider the family of 
matrices with columns indexed by S whose row space equals L. The matroid is presented 
by the linear dependencies among the columns of any such matrix. When L is viewed as a 
linear subspace of functions from S to K, the matroid is the "function space geometry (or 
chain-group geometry) G(S, L)" discussed in [63, §1.1.C] . As such, each e £ S corresponds 
to the linear functional given by evaluation of / £ L on e. Therefore these functionals, as 
a finite subset of the dual space of L, comprise a vector representation of the matroid. 

The members of L present an oriented matroid 7V(N) by a set of covectors C. Each 
covector is the function I : S —> {+, — , 0} determined by some / £ L by /(e) = sign(/(e)) 
for all e £ S. In other words, if L is viewed as the row space of a matrix, then each 
I £ L is the signature (i.e., the sequence of +, — or 0s) that indicates the signs in one row 
/ £ L. The signed cocircuits are the covectors with minimal non-empty support. See [1] 
for an exposition of oriented matroids that begins with linear subspace presentations 
including the cycle and coboundary (or cocycle) spaces of graphs. Our topic utilizes 
exterior algebra and the chirotope given by the signs of the Pliicker coordinates of N to 
present oriented matroids. Theorem 12.11 states the needed details. Deeper discussions 
appear in [4, especially §2.4 on stratifications of the Grassmann variety and chap. 8 on 
realizations]. 

The following theorem summarizes the facts we will need. It characterizes those ele- 
ments of S(V) (as fully decomposable) that determine linear subspaces and present the 
.fC-realizable matroids with ground set S. We will call such elements extensors. 

Theorem 2.1. Given non-zero element N £ £{V) where V = KS, the following condi- 
tions are equivalent: 

1. There exist r linearly independent vectors Vi £ V for which 

N = vi ■ • • v r . 

(When r = 0, the condition is N = al for some a £ K , a ^ 0.) 
L is the subspace spanned by {v\, . . . , v r }. 
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2. There exists r, < r < \S\, such that the only non-zero coefficients N[A] in 

N = ^ NL4]A 

ACS 

satisfy \A\ = r, and the function NL4] /rom sequences A to K is alternating and 
satisfies the Grassmann-Pliicker relationships: 

For all length r sequences A = a\ ■ ■ ■ a r and B — b\ ■ • ■ b r over S , 

r 

N[A]N[B] = ^N[b i a 2 ---a r ]N[b 1 ---aibi---b r ]. 

i=l 

(Here, afci means bi within sequence B is replaced by ai.) 

3. There exists a rank r matrix N with r rows and with columns indexed by S for which 
the coefficients o/N with \A\ = r satisfy 

N[A] = detN(A) 

where N(A) is the submatrix of N with columns A, and the other coefficients are 0. 
(For r = the condition is N[0] ^ and all other N[A] = O.j 

L zs i/ie subspace spanned by the rows of N. 

J/r ^ then N and the v\,...v r in N = Vi • • • v r can be chosen so row i of N holds 
the coefficients for writing Vi as a linear combination of vectors from basis S . 

Definition 2.2 (Extensor). When the conditions in Theorem 12.11 are true about N, 
we say N is decomposable 2 , and the coefficients denoted by N[A] are called the Pliicker 
coordinates of N. A decomposable element of an exterior algebra is called an extensor. 
The integer r is called its rank and is denoted by pN. 

Remark: We do not define Pliicker coordinates to equivalence classes of homogeneous 
coordinates, so extensors that differ by a non-zero scalar multiple have different Pliicker 
coordinates, even though they represent the the same subspace. 

Remark: NLA] is defined as for all A with \A\ ^ pN. 

The fact that each rank r extensor is the exterior product of r vectors, and V1V2 = 
— V2V1 for vectors implies that extensor multiplication satisfies the following anticom- 
mutative law: 

NiN 2 = (-l) pNlpN2 N 2 N!. (4) 

Theorem 2.3. Suppose r, subspace L, matrix N , extensor N(S) and its Pliicker coordi- 
nates are as described in Theorem \2.1\ 

2 We follow [31] and other authors who omit the qualifier "fully" in this context. 



13 



1. There exists a rank \S\—r matrix A r± with columns indexed by S and with \S\ — r 
rows for which coefficients with \A\—r satisfy 

N[A] = det A^(A)e(A4) 

where A — S \ A expressed in an arbitrary sequence and e is some non-zero alter- 
nating sign function of sequences over S . 

2. L consists of all x G KS that satisfy the equations 

N L x = 0. 

In other words, L ( the row space of N) and the row space of A r± are orthogonal 
complements in KS. 

Proof. See [31, chap. 7]. An elementary proof of theorem 12.31 related to theorem 12.11 
appears in [30, VII. 3 Theorem I]. The equivalence of the Grassmann-Pliicker relationships 
to the other conditions is proved in [30, VII. 6 Theorem II]. □ 

An extensor N with ground set T is the finite set T paired with an extensor in 
8{KT). We use the notation N = N(T) and T = S(N) to indicate that N has ground 
set T. 

We need the ground set for the same reason that a ground set is necessary to define 
the dual of a matroid with coloops. Independent sets, or the collection of sequences B for 
which chirotope x{B) = + is not sufficient because a loop doesn't appear in any of these 
objects. Furthermore, in our calculations and proofs we find it very helpful to combine 
the signs of sequences (i.e., permutations of subsets) from several different sets within one 
analysis without having to relabel any elements. Identities like the Tutte equations relate 
function values for objects with different ground sets. It is not sufficient for an extensor 
or chirotope to be defined up to sign for certain identities to be valid (not just up to sign). 
This validity may facilitate the use of the identities in computer programs. 

Theorem 2.4. Given N(jS'), let N be a matrix satisfying Theorem \2. li 

1. The collection of those B C S for which N[JB] ^ is the collection of bases of a 
matroid with ground set S. 

The same matroid is presented by the independent sets of columns of N. 

2. The function \ of sequences over S into {+1, —1,0} for which x(B) ^ s the sign of 
N[£?] is the chirotope function of the oriented matroid M denoted by M = Af(S) = 
JV(N). 

The covectors of M are presented by the signatures of N 's row space elements; 
the signed circuits (i.e., oriented matroid "vectors" with minimal support sets) are 
presented by the signatures of the minimal linear dependencies among the columns 
ofN. 
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3. If ~N[B] £ {0,±1} for all B, then 7V(N) is the unimodular (or regular) oriented 
matroid whose chirotope function satisfies x{B) = N[S]. Furthermore, every uni- 
modular oriented matroid can be presented by such an N. 

Proof. See [4] or [1, chap. 5]. Details pertaining to the unimodular matroids including 
several characterizations are given in [64, Theorem 3.1.1, p. 41]. □ 

Definition 2.5. If N(S) ^ and e £ S then 

• e is called a loop if N[JB] =0 for all B with e £ B, and 

• e is called a coloop if every B such that N[.B] ^ satisfies e E B. 

Remark: e is therefore a loop or coloop in N if and only if it is a loop or coloop 
respectively in the matroid presented by N. 

Definition 2.6. Each N(S') 7^ defines the function pn on subsets ACS where Pn(A) 
is the rank of A in the matroid presented by N(S'). 

Theorem 2.7. Given N(S) ^ 0, e £ S and S' = S \ e: 

1. The Pliicker coordinate function for N restricted to sequences B C S' is the Pliicker 
coordinate function for an extensor denoted by (N\e)(S"). This operation N — > N\e 
is called deletion of e. 

(N \ e) 7^ if and only if e is not a coloop in the matroid presented by N. In this 
case, the unoriented or oriented matroid minor M \ e is presented by N \ e and 
p(N\e) = pN. 

2. The function defined by N[i?e] for sequences B C S' is the Pliicker coordinate 
function for an extensor denoted by (N/e)(5"). This operation N — > N/e is called 
contraction of e. 

(N/e) ^ if and only if e is not a loop in the matroid presented by~N. In this case, 
the unoriented or oriented matroid minor N/e is presented by N/e and p(N/e) = 
pN-1. 

Remarks: If N = then N/e = N\e = 0. The zero extensor does not present 
any matroid. All rank (empty S or loops only) matroids have only one basis 0; they are 
presented by the non-zero extensors al(S) of rank 0. 

Proof. Let iV be a matrix representing N in Theorem 12.11 

(N \ e)(5") is the extensor known from Theorem 12.11 when the column labeled by e 
is deleted from N. Note that if e is a coloop then this reduces the rank of N and so 
(N\e)(S') = 0. 

If e is a loop in N then (N/e) = 0. Otherwise, the N[Se] are fixed non-zero multiples 
of the Pliicker coordinates from a matrix obtained from the N by row operations to make 
all but one entry in column e zero and then deleting the row and column with that 
non-zero entry. 

See [4, §3.5] for oriented matroid minors and other structures in terms of chirotopes. 

□ 
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Theorem 2.8. Given N(S) and e G S, 

N(5) = (N/e)e+(N\e)l({e}) 

The multiplication by l({e}) makes the ground set of the second term be S instead 
of S \ e. It will be omitted in contexts where the ground set is clear. 

Proof. Let B C S. We prove that each Pliicker coordinate N[B] equals the sum of the 
corresponding Pliicker coordinates of the extensors on the right. 

If e G B we can write B = B'e. N[B'e] = (N/e)[5'] = (N/e)e[S'e], and (N\e)[S] = 0. 

lie&B then (N/e)e[S] = and N[B] = (N \ e)[B\. □ 

It is convenient to let N/A denote N/a^/ ■ • • ja\ where A — a\ ■ ■ ■ a^, and similarly for 
N \ A. It follows that N/A[X] = N[XA) for all X. We note that for a G 6 fc , 

N/A a = e(a)N/A, but 
N \ A a = N \ A. U 

2.3 Ground Set Orientation and Duality 

Definition 2.9 (Ground set orientation). An orientation of the ground set e is an 
alternating function into {+1,-1,0} of sequences of ground set elements that is non-zero 
on sequences of distinct elements, and which satisfies e(0) = 1. 

One family of ground set orientations is derived from fixed linear orders on all ground 
set elements using the rule that e(X) = (— l) v where v is the number of inversions in 
X (where an inversion is with % < j and x% > xj). A permutation a G & n of 

{1, . . . , n} is always considered a sequence o\o~i . . . o~ n of natural numbers with ground set 
orientation derived from their usual ordering. Hence, e(o~) is the usual sign of permutation 
a. However, ground set orientations of matroid elements or graph edges will not be 
assumed to derive from a linear order. 

Since permutations a G & n and sequences of ground set elements will not be confused, 
we use the same symbol e for permutation sign and ground set orientation. 

Given a length n sequence X — x± . . . x n and a G & n , let X a denote x ai . . . x an . The 
following routine facts will be used in our proofs: Of course, F is alternating means F(X a ) 
= e(cr) F(X) for all sequences X and a G &\x\- 

Lemma 2.10. Suppose e\ and e 2 are arbitrary alternating functions of sequences. 

1. If n = \X\ = \Y\, a G & n , and A, X, B , C,Y,D are sequences then 

e 1 (AXB)e 2 (CYD) = e(a)e 1 (AX a B)e 2 (CY D) 

= e(a)e 1 (AXB)e 2 (CY a D) 
= e x {AX a B)e 2 {CY a D) 

2. e i {XY) = {-l)\ x W Y k i {YX). 
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With a ground set orientation e in hand, we define: 

Definition 2.11 (Canonical Dual). Given N(S), N X [X] = N ±S LY] = N[X]e(X X), 
where X is any sequence of the distinct elements in 5* \ X. 

The symbol _L e will be abbreviated by _L when e is irrelevant or doesn't require em- 
phasis. 

Remark: Each ground set S determines with e a sign choice from among the two 
that both provide a presentation of the oriented matroid dual. 

The demonstration in [4, end of §3.6] of oriented matroid chirotope dualization has a 
similar formula, whose right hand side is independent of an arbitrarily chosen sequence. 
It follows that our extensor dualization corresponds to the oriented matroid dualization 
of the oriented matroid presented by N. Theorem 12.31 justifies this for realizable ma- 
troids. (Dualization is also the Hodge star operator [30] when S is identified with the 
corresponding basis of the dual space. Also see ^7. II ^ 

2.4 Identities 

Our main proof uses some identities on extensors that correspond to well-known relation- 
ships among matroid operations. These identities involve extensors with ground sets for 
which a ground set orientation is used to define dualization. The union of disjoint sets is 
denoted by U . 

Theorem 2.12. 

(Ni + N 2 )- L = Ni- + N£ 

(aN) 1 = aN 1 1 J 

N LX (S) = (_i)pN(|shpN) N(5) = (-l)" NpN± N(S) (7) 
Given N(S'), and sequences X C S and S' = S \ X , 

(TSt\X) ± = e(S')e(S'X) (N x /X) (8) 

(N/X) ± = e(^ / )e(5 / X)(-l) |x|(|5| -" N) (N x \ X) (9) 

Given Nj(S'j) with SiflS^ = 7 the extensor product NiN2(Si US2) presents the (oriented) 
matroid direct sum and 

(N1N2) 1 = e(S l )e(S2)e(S 1 S 2 )(-lY^^ (10) 

Proof. Linearity (jUJ) is immediate from the Pliicker coordinate definition. It will be re- 
peatedly used with a = ±1 below. 
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To prove Theorem 12.121 (J7J) , write 

N ±J -[A] = (N^iA] = (N^A] e(AA), = N[A] e(A A)e(AA) = N[A] e(AA)e(AA) 

where in the last equation we chose the sequence order A = A. Therefore the sign 
correction is (— 1)' A ' For non-zero coordinates this is (— l) pN (l 5 ! - ^) = (_i)p n p n . 
To prove Theorem 12.121 (JSJ) , write 

(N \ X^iA] = (N \ X)[A] e(AA) = N[A] e(AA), 

where A — S'\A. But in 

(N ± /X)[A] = N^AX] = N[AX] e(AXAX) 

the elements in sequence AX are S' \ A, the same as in the sequence symbolized by A in 
the previous equation. We can therefore choose AX = A and write 

(N"7J0[A] = N[A] e(AAX). 

Combining the two sign corrections gives e(AA)e(AAX). That equals e(S')e(S'X) for all 
reorderings 5" of AA. 

We can get © from Q and ©. Specifically (N/X) 1 = (L ± /X) ± (-1)^ nX with 
L = N" 1 . This equals 

(L \ X) x± e(S')e(S'X)(-l) pNpN± 

= (N X \X) (-l)"( N± \ x ^ N / x )e( 1 S / )e(.S / X)(-l)^ N± 

As usual, we can restrict attention to non-zero coordinates. Let pN = r, \S\ = s, and 
\X\ = x so p(N/X) = r — x and p(N ± \ X) = pN -1 = s — r. The sign correction is 
therefore 

€ (S')e(S'X)(-lY s - r){r - x)+r{s - r) = e(S')e(S'X)(-iy s ~ r){2r - x) = e(S')e(S'X)(-l) {s - r)x . 

To prove Theorem 12.121 (fTUj) . take A;lA 2 = Ai A 2 with each Aj = Si \ Aj in 

(N 1 N 2 ) ± [A 1 A 2 ] = 
(N!N 2 )p^ e(CTA!A 2 ) 
= N 1 [A7]N 2 [A^]e(A7A 2 ~A 1 A 2 ) 
= N^Ax] e(ATA0 N 2 L [A 2 ] e(^A 2 ) e^A^A^ 
= N^[Ai] N 2 L [A 2 ] e(A 1 "A 1 )e(A;A 2 )e(A 1 "(A 1 A;)A 2 )(-l) |Al| l^l 

= N^AJ N^A,] 6(5 1 )6(^ 2 )6(5 1 5 2 )(-1)" N ^ N2 

where in the last equation, we applied permutations a and r, each twice, for which 
(AiAi)^ = Si and (A 2 A 2 ) T = S* 2 . We then substituted the correct ranks for cases where 
the coordinate is not 0. □ 
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3 An Extensor Tutte Function 



Recall that a ported extensor or matroid is one whose ground set has a distinguished 
subset of port elements. 

Given a ported extensor N(P, P) (the notation means P is the set of ports and the 
ground set is P UP), we will define a parametrized extensor M£(N) using extensor 
operations. We will illustrate its construction with extensors and equivalent matrices. We 
then write and prove parametrized identities satisfied by the function N(P, E) — > M£(N) 
which are analogous to the ported Tutte equations. Our identities however apply to 
extensors rather than to commutative ring values. The identities include sign- correction 
factors that depend on the particular ground set orientation e that was used to define 



The definition of M^N) below applies to all extensors N over K{E U P). The main 
result therefore belongs to exterior algebra. Section 0] shows how defines the extensor 
valued function on the minor closed class of ported unimodular oriented matroids Af(P, E) 
by M^A/") = M£(±N) where M is presented by either unimodular extensor ±N. 

Each linear map on V can be extended to a unique exterior algebra map on the exterior 
algebra S(V) [31, Theorem 7.1]. Given ground set P U E, let P v and P L be two disjoint 
copies of P, both also disjoint from E. For each p G P, let p v G P v and p L G P L be 
the corresponding elements in the respective copies. We define the following maps from 
K{P U E) to K[r e , g e ](P v U P L UP) and extend them to the exterior algebra. If necessary, 
the field K is extended with the parameters g e ,r e , e G E. 



In terms of matrices, i g signifies multiplying column labeled e by g e for each e G E 
and renaming column p by p b for each peP. Likewise, v r signifies multiplying column e 
by r e and renaming column p by p v . 

The parameter subscript symbols in i g and v r will sometimes be omitted for brevity. 
For subset Q C P, Q v denotes {q v : q G Q} C P v and Q L denotes {q L : q G Q} C P L . 
Recall that set symbols denote sequences. The sequences of Q v and Q L correspond to the 
sequence of Q. 

Definition 3.1. Given a ported extensor N(P, P), a ground set orientation e and dual 
operator _L e , parameters g e and r e for each e G P, and e-preserving functions v r and t g 
defined above, let 



Hence, M E (N) is defined as a ported extensor M(N) = M(N)(P t U P v , E) contracted 
by the sequence of non-port elements P. Therefore, by Theorem 12.71 it is an extensor. 
Each pair of sequences / C P, V C P with \I\ + |V^| = |P| specifies the Plucker coordinate 
of M£(N) with index I L V V and value 



M B (N). 



v r (e) = r e e for e G P and v r (p) = p v for p G P. 
i 9 (e) = (7 e e for e G P and £ 9 (p) = p t for p E P. 



(11) 



M(N) = t ff (N) w r (N ±e ) and 
Mb(N) = M(N)/P 



(12) 



M S (N)[PK] = (t(N) ^(N^))^^]. 



(13) 
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Figure 1: Graph defining the graphic oriented matroid M 
Proposition 3.2. For a E K, M E (aN) = a 2 M E (N). 

Proof. M(aN) = a 2 M(N) is immediate from the definition. Contraction "M./E is linear 
in M. □ 

We can express ()12)) in matrix terms. Let N be some full row rank matrix with 
columns indexed by P U-E that presents N(P k)E). Similarly, let N denote a matrix 
presentation of N -1 . 

Example. We show one totally unimodular matrix representation N of the ported 
graphic oriented matroid with P = {pi,P2,P3} and E = {ei, e 2 , e 3 , e 4 } for the graph in 
figure n The rows code 3 oriented cutsets which determine a basis for the 1-coboundary 
(or cocycle) space. We also express N by the exterior product of the vectors given by the 
rows of this matrix. 



N 



Pi P2 P3 

-1 +1 

+1 -1 

-1 -1 +1 



ei e 2 

+ 1 +1 
-1 
+ 1 



e 3 









e 4 


hi 
+1 



(-Pi + P3 + ei + e 2 )- 
N = (P2 - P3 - ei + e 3 )- 

(-pi - p 2 + p 3 + ei + e 4 ) 

Next, we write one totally unimodular matrix A r_L for the canonical dual. We have 
checked that the sign satisfies Definition 12. Ill with e chosen so e(piP2P3Cie2e 3 e4) = 1 by 
verifying N- L [eie 2 e 3 e 4 ] = N[pip 2 p 3 ]e(pip 2 P3eie2e 3 e 4 ). 



Pi 


P2 


P3 


ei 


e 2 


e 3 


ei 








+1 


-1 











+1 


+1 


+1 











+ 1 





+1 


+1 





-1 








+1 





+1 








+1 






iV J 



Continuing the general discussion, let G and R be the diagonal matrices of the g e and 
r e . The matrix 

" N(P) N(E)G 
N ± (P) N ± (E)R 



M(N) 



(14) 
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has order (p + e) x (2p + e), columns indexed by sequence P L P V E it and presents M(N). 
Example continued. We abbreviate labels p L \ and p v i by i\ and f 1; etc. 



M(N) 



H 


1 2 




Vi 


v 2 


^3 


ei 


e 2 


e 3 


e 4 


-1 





+1 











9i 


92 











+1 


-1 











-9i 





#3 





-1 


-1 


+ 1 











9i 

























+1 


-n 




















+1 


+1 


+1 











''1 














+1 


+1 





-r 2 

















+1 





+1 








r 3 






Generally, M(N)(P t U P v U E) is the exterior product of the vectors in K(P L UP„UE) 
corresponding to the rows of this matrix. M^(N)(P t U P v ) appears in the expression 

M(N) = (M B (N))E+--- 

where the initial term is the only one with factor E. 

Example continued. We calculate M^(N) by doing ring operations on rows to 
eliminate all but one non-zero entry in each E column in M(N). The result is that 

gig 2 g3g 4 rir 2 r 3 r 4 M(N) 

is equal to the following extensor in matrix form: 





1 2 




Vl 


v 2 


^3 


ei 


e 2 


e 3 


e 4 


—T\T2 










air i 


9\T<1 + 92T\ 

















rir 3 


-nr 3 


-gsn 





-9in - gsri 














—r\r 4 


—r 1 r 4 


r±r 4 


-94?! 


-94,n 


gir 4 - g 4 r x 





























9i 


-gin 




















9aXi 


9aTi 


g^Ti 











g 4 r x r 4 














9ir\ 


g 2 r\ 





-g 2 r 1 r 2 

















9sri 





gi.fi 








9arir 3 






After some cancellation, we can read off the answer from the 3x6 upper left submatrix, 
which is a matrix presentation of the extensor rfM^(N): 



■'1 


'2 


«3 




t'2 


V 3 


r±r 2 





T\T 2 





g2n 


gir 2 + g 2 r\ 





rir 3 


-nr 3 


-g^n 





-gir 3 - gzn 


T\T 4 


—r 1 r 4 


r x r 4 


-94Ti 


-94T1 


gir 4 - g 4 r l 



Remark: Each Pliicker coordinate of M^(N) is a homogeneous polynomial of degree 
\E\ in the g e , r e . However, this example demonstrates that there sometimes doesn't exist 
a matrix expression for M^(N) all of whose entries are polynomials. The reader can 
verify that each order 3 minor of the above matrix is a multiple of r\. 
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Graph Minor Extensor Minor Term in M^(N)[viiiV3] 




N/{ei,e 4 }|P -gmr 2 r 3 



N/{e 2 ,e 3 }|P +9293^4, 



Figure 2: The two graph and extensor minors with corresponding terms in M_e(N) [^liii^]. 

Example continued: Here are examples of Pliicker coordinates, which can be cal- 
culated from the above matrix as order 3 minors divided by r\. See §6.31 

M E (N) [v x v 2 v 3 ] = g x g 2 g 3 r A + 9x929^ + 9\9z9^2 + 9i9z9^r x 
M E (N)[iif 2 w 3 ] = {gin + gzri)(g 2 rA + g±r 2 ) 

M E Q$)[viiiV 3 ] = ~gmr 2 r 3 + g 2 g 3 rir 4 

Observe Me(N) [vi^^] is the basis enumerator for J\f(N) \P. The graph and extensor 
minors corresponding to the terms of M^(N)[f ii\V 3 ] are shown in figure El 

3.1 Main Result 

Theorem 3.3. The parametrized extensor-valued function M^(N)(P„ UPJ of ported ex- 
tensor N = N(P, E) has the following properties: 

1. Given Ni(Pi, Pi) and N 2 (P 2 , E 2 ) with E = E x U P 2 and P = P x U P 2; 

M E (N! N 2 )(P,P)=e(P 1 P 2 P)e(P 1 P 1 )e(P 2 P 2 )M Sl (N 1 ) M S2 (N 2 ). (15) 

2. If e C E and E' = E\e then 

M E (N) = e(PE)e(PE') (<7 e M fi ,(N/e) + r e M B ,(N \ e)). (16) 

3. Let E = 0. The Pliicker coordinates o/M (N)(P t UP„) satisfy 

M (N)[J t K] = M[I L V V ] = e(V V) N[J]N[F]. 
for all I CP and V CP. 

4. M E (0) = 0. 
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Properties 1. and 2. in terms of Pliicker coordinates are: 

M E (N 1 N 2 )[I u V lv I 2I V 2V ] = 
e{P 1 P 2 E)e{P 1 E 1 )e{P2E 2 ) M. E ^i)[hiV lv \ M Ea (N 2 ) [I 2I V 2V ] . (17) 

M E (N)[iy v ] = 

e(PE)e(PE') [{g e M. E i (N/ e) [I b V v ] + r e M. E t (N \ e) [I t K] ) ■ (18) 

Remarks: 

1. Property 2. implies that every linear combination gJs/l E i{N / e) + r e M^(N \ e) is 
decomposable, i.e., an extensor. 

2. Proposition HO with a = ±1 implies M^NiNa) = M jB (N 2 Ni). We can also verify 
this from the right hand side of property 1 using Lemma (J2.10I 3) and noting the 
rank of extensor Mg^N,) is |Pj|. 

3. If N ^ 0, one but not both of N/e and N \ e will be the extensor if and only if e is 
a loop or a coloop in the matroid of N. If N' = then N /J - = and M^(N / ) = 0. 
We therefore write property 2 without restricting e to a non-separator. 

4. Property 1 except for signs is immediate from direct sum of subspaces and their 
corresponding extensors. Property 2 except for the signs follows immediately from 
the fact that minor [I t V^,i?] of matrix ()14|) equals a linear combination with coeffi- 
cients g e and r e because the column e belongs to this minor no matter which e G E 
is specified for the identity. 

Proof. From the definition, M(NiN 2 ) = t(NiN 2 ) v((NiN 2 )- L ) which equals 
e(S 1 )e(S 2 )e(S 1 S 2 )(-iy^^ t(NiN 2 ) v(N^) 

by Theorem l2~T2¥fll . t(NiN 2 ) u(N^N£) = t(Ni)t(N 2 )v(Nf- )f(N 2 L ) which equals (by 
O) 

(_1)PN^ P N 2 ^ Ni) V ( N ±^ v ( N ±y 

Therefore M(NiN 2 ) = 

e(S 1 )e(S 2 )e(S 1 S 2 ) t(N0 u(N^) t(N 2 ) u(N^). (19) 

Therefore, M(NiN 2 )/£i£ 2 = 

e{S 1 )e(S 2 )e(S 1 S 2 ) ((M(N 1 )/ J E 1 ) E x (M(N 2 )/£ 2 ) E a )/^£? 2 

because the only Pliicker coordinates of the form M(Nj)[JQ] for i = 1 or 2 that contribute 
to (|19|) when it is contracted by E\E 2 satisfy E{ C Xj. The anticommutativity law Pjl 
then implies that M(NiN 2 )/£i£ 2 = 

e(5 1 )e(5 2 )e( 1 S 1 5 2 )(-l)l £;i ll P2 l (M^Nx) M Sa (N 2 ) E x E 2 )/£ 1 £ 2 
= e(5 1 )e(5 2 )e(5 11 S 2 )(-l)l £;i ll P2 l M^Nx) M B2 (N 2 ). 
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Since the sequence orders of the Si are arbitrary, let Si = PiEi for i — 1 and 2. According 
to equation ©, M^N^) = M(N 1 N 2 )/ J B = e(a)M(N 1 ~N 2 )/E 1 E 2 where {E X E 2 )„ = E. 
Therefore, M E (NiN 2 ) = ±M Bl (Ni) M£ 2 (N 2 ) with the sign equal to 

e(a)e(P 1 E 1 )e(P 2 E 2 )e(P 1 E 1 P 2 E 2 )(-l)\ E ^ 
= e(a)e(P 1 E 1 )e(P2E 2 )e(P 1 P 2 E 1 E 2 ) 
= e 2 (a)e(P 1 E 1 )e(P 2 E 2 )e(P 1 P 2 (E 1 E 2 ) a ) 
= e(P 1 E 1 )e(P 2 E 2 )e(P 1 P 2 E), 

which proves property 1 of the theorem. 

Now for property 2. Let us apply Theorem 12.81 to N and N -1 , and apply i and v 
respectively. 

N = (N/e)e+ (N\e)l(e). 

N x = (N ± /e)e + (N ± \e)l(e). 

l(N) =t((N/e)e) + t ((N\e)l(e)) 
=i(N/e)^e + t(N\e)l(e). 

V (N X ) =v((N ± /e)e) +v((N ± \ e)l(e)) 
=u(N ± /e)r e e + v(N ± \ e)l(e). 

The exterior product of t(N) and t^N -1 ) is therefore 

g e t(N/e) e ^(N x \ e) + r e i(N \ e) v^/e) e 
+# e r e t(N/e) e u(N J -/e) e 
+i(N\e) ^(N^e). 

M^(N) is the result of contracting the above by £7. The third term above is because e 
is a repeated factor. The last term will vanish when contracted by E because none of its 
non-zero Plucker coordinates have an index that contains e. So we will omit them in the 
following. By Theorem 12 . 1 2t| §j) 

^(N 1 - \ e) = v((N/e) ± ) e(5 / )e(5 / e)(-l) 1 -^l^ N ) 

and 

u(N-7e) = v((N\ e) X ) e(S")e(S"e)- 

Notice that p(N/e x ) = \S'\-(pN-l) = \S\-pN. So e ^(N/e)^ = ^(N/e)^ e(-l)l 5 l^ N . 
Therefore when the above substitutions are made we get 

M E (N)e(S')e(S'e) = (^(N/e) ^((N/e) 1 ) + r e .(N \ e) u((N \ e) 1 )) e/£. 

Lemma f2. 101 used with a such that E a = E'e shows that the right hand side is 

e(a)((...)e)/E'e = e(E)e(E'e)((---)e)/E'e. 
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So the right hand side is 

e(E)e(E'e)(g e ME'(N/e) + r e M B ,(N \ e)). 

Since the sequence order of S' is arbitrary, we can choose S' = PE' . The sign correction 
is then 

e(S')e(S'e)e(E)e(E'e) = 
e(PE')e(PE'e)e(E)e(E'e). 

Applying the permutation r for which (E'e) T = E to the two appearances of subsequence 
E'e does not change this expression's value. Hence the sign correction is 

e(PE')e(PE)e(E)e(E) = e(PE')e(PE) 

and property 2 of the theorem is verified. 

The definition of Mg immediately gives property 4, and, together with the definition 
of extensor dual, gives property 3. □ 

Corollary 3.4. The set o/Mg(N) obtained as the g e , r e range over M for each e G E 
represents the points in a projective subspace of a Grassmannian (which consists of all the 
linear subspaces overM>(P L UP„) with dimension \P\). 

Proof. Induction: Use Theorem 13.31 property 2 for when \E\ > and property 3 for the 
the basis. □ 

Proposition 3.5. Given N = N(P, E), and sequences I C P, V C P, and V = P \ V , 
e(V V)e(PE)M E (N)[I L V v ] = e(P) £ N[IA]N[VA]g A rj. 

ACE 

Remark: The only non-zero terms in this sum are those for which both A U / and 
AUK are bases in the matroid of N. 

Proof. Recalling that E symbolizes a sequence e\ ■ ■ -e n , let E^ = e^j+i so E\ — E and 
E n+ i = 0. When property (2) of Theorem 13.31 is applied successively for e := e^, E := Ei 
and E' := E i+ \ for i — 1, 2, . . . , n and the products are expanded, the result is a sum of 
2 n terms, one for each subset A C E. For each fixed i, 1 < i < n, the instances of the 
symbols E and E' — -E 1 \ e, within all applications of property (2) each denote the same 
sequences. Therefore, sign cancellation occurs and we can write 

Mb(N) [hVv] = e(PE) e(P)g A r A M d (N /A \ A)[I b V v \. (20) 

ACE 

Property (3) combined with the definitions of extensor deletion, contraction and dualiza- 
tion demonstrate that within each term 

M (N/A \ A)[I L V V ] = e(VV)N[IA]N[VA], 

and the conclusion follows. □ 
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Corollary 3.6. If N 7^ and the parameters are generic or are all positive, then 
Mjj(N) ^ 0. 

Proof. N has some basis, i.e., B C P U£ for which N[P] 7^ 0. (N.B. B = is possible.) 
Take J = 5 n P and V = P \ I. Proposition E3] indicates M E (N) [I L V V ] ^ with A = 
B\I. □ 

Corollary 3.7. Pliicker coordinate M£;(N)[J t V / u ] is a homogeneous polynomial in the g e , 
r e whose terms are square-free and have degree pN — |/| in the g e and degree \E\ — pN + \I\ 
= \E\ + |P| - pN - \V\ = pN- 1 - \V\ in the r e . 

Proof. Immediate from Proposition 13.51 matroid duality and the fact N[X] 7^ only if 
\X\ = pN. □ 

Corollary 3.8. 

e(PE)M E (N) = e(P) ^ M (N/A|P)^. (21) 

A C E : p N A = \A\, 
pN - p(N/A|P) - p N A = 

Proof. The definition of extensor deletion indicates N/A|P is an alternative notation for 
N/A \ A when A <0 E, A means E \ A and the ground set of N is P U E. The conditions 
stated in Corollary 13.71 allow us to restrict the sum as indicated. Hence formula (|2T?j) for 
each Pliicker coordinate is equivalent to the given expression for the extensor. □ 

The following definition and consequence of Proposition 13.51 clarify some of the sign 
behavior resulting from the definitions. 

Definition 3.9. A function F = F e (X) whose value might depend on the ground set 
orientation e and on the sequence X is said to be 

1. alternating in X if F e (X a ) = e(a)F e (X), for all a G &\x\] and 

2. alternating in e if F~ e (X) = -F e (X). 
Corollary 3.10. Let Q C P L U P v with \Q\ = \P\. 

1. M^(±N)[Q] is constant under sign change of ±N, and is alternating in E, e and 
Q. 

2. e(PP)M|;(±N)[Q] is constant under sign change of ±N and under changes or 
reorderings of e or E; it is alternating in P and in Q. 

3. e(P2£)M^(±N)[P t ] enumerates the bases o/jV(N/P), assuming P is independent 
in the matroid jV(N) ; by 



e(PE)M|(±N)[Pj = ]T^N 2 [PP], 



BCE 
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4- and e(PE)~M. e E (±'N)[P v ] enumerates the bases ofAf(N\P), assuming P is coinde- 
pendent in 7V(N) ; by 

e(PE)M%(±N)[P v ] = 9bTb^ 2 [B]. 

BCE 

Remark: Properties 3. and 4. express the Matrix Tree Theorem. 

4 Corank- Nullity Polynomials 

The well-known corank-nullity (or rank) polynomial is easily generalized to ported ma- 
troids. We did [12] this with a definition that differs from Las Vergnas' big Tutte polyno- 
mial [37] only in notation and in our applications. We will now generalize it by including 
parameters and modify it for oriented matroids by reinterpreting the symbols for minors. 
In the definition below the symbol J\f/A\P represents the oriented minor of the oriented 
matroid M obtained by contracting A and restricting to P. Minors that are different 
orientations of the same unoriented matroid are deemed different objects. 

Definition 4.1 (Parametrized and Ported Corank-Nullity Polynomial). 

R(Af(P,E)) = £ [N/A\P] g A r^u pK ~ p ^' A \ p ^ pA v^~ pA . 

ACE 

In this formula, the bracketed oriented matroid [J\f/A\P] = [Afn © ... © Af ic ] denotes 
the (commutative) product of the variables [J\fa] , . . . , [N"i C ] , where each variable signifies 
a connected component of Af/A\P. If P = then [J\f/A\P] = [0] = 1; so R reduces to 
the corank-nullity polynomial, parametrized. 

The formula therefore defines a polynomial in parameters g e , r e , whose (other) variables 
are u, v together with a distinct variable for every connected component of every minor 
of M obtained by contracting some subset A C E and deleting A — E \ A. The latter 
variables only occur in monomials that signify direct sums of one or more minors. 

It is readily verified that R(J\f(P, E)) satisfies the ported Tutte equations below. The 
details published in [12] can be immediately adapted to the changes we described above. 
We state these results without proof: 

Proposition 4.2. 1. If e G E is neither a port nor a loop nor a coloop in J\f(P, E), 

R(M(P,E)) = g e R(N/e)+r e R(N\e). (22) 

2. R(M l ®M 2 ) = i?(M)^(Ar 2 ) = i?(Ar 2 )i?(M). 

3. R(N\{e)) = g e + r e u and -R(Ao(e)) = r e + g e v, for the coloop and loop matroids 
M(e) and W (e) on E = {e}, P = 0. 

I R(Af(P, 0)) = [A/] (i.e. when E = 0.; 
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Let us take M to be the oriented matroid presented by extensor N. The reader can 
now verify that the expansion in corollary 13.81 for M#(N) is obtained from R(Af(P,E)) 
by substituting u = 0, v = and the extensor e(P)e(PE)'M.^(N / A\P) for the monomial 
[J\f/A\P] in the term with factor g^rj, for each A. Note that \A\ — pA = and pM — 
p[M / A\P] — pA = pN ' — p(P U A) = imply that A is independent and P U A is spanning 
in M . Therefore N/A\P ^ for those terms where the exponents of u and v are both 
zero. 

With arbitrary N, the substitution of e(PP)e(P)M (N/,4 l |P) for monomial [M/A^P] 
= [J\f/Ai\P] in R(J\f(P,E)) is not well-defined. The reason is that the same oriented 
matroid M /Aj | P might be represented by different extensors all with the form N/Aj|P, for 
various Aj, ^ Aj . They may differ by representing different subspaces (in the same oriented 
matroid stratification [4, §2.4] layer) of KP. Therefore, different values M£(N / Ai\P) 
must be substituted in [M / Ai\P]g a^^- with different Ai even though these [J\f/Ai\P] all 
denote the same oriented matroid. 

The one general situation where R(J\f) with u — v — determines M£(N) is when 
N(P, E) is a unimodular extensor, i.e., one that represents the unimodular oriented ma- 
troid M . [64, Theorem 3.1.1, p. 41] provides this among other equivalent characteriza- 
tions of unimodular (also called regular) matroids. One of these characterizations is that 
bracket values from {+1,-1,0} may be assigned so the Grassmann-Pliicker relationships 
hold over Q. 

Definition 4.3. The extensor-valued function Af(P,E) — > M^A/") is defined on the 
minor-closed class of ported unimodular oriented matroids by 

M £ (AO = Mg(N) 

where ±N(P, E) are the two unimodular presentations of the ported oriented matroid 
Af(P,E). 

The simplest case to demonstrate that the unimodular matroids must be oriented in 
order for monomial substitution in P(A/"(P, P)) to produce M^A/") is the two orientations 
Afi, A/2 of the 2-circuit matroid on two ports. Here, E = and M0(A/i) 7^ Ma (A/2). 

We conclude: 

Theorem 4.4. For unimodular oriented matroid M = Af(P,E), M^A/") is the result 
of evaluating R(Af) (in the exterior algebra) after the substitutions u = 0, v = and 
extensor e(P)e(P£)Mg(A/i(P)) for each monomial [A/"j(P)] (which symbolizes an oriented 
unimodular matroid with ground set P). 

Proof. Immediate from the above remarks and Corollary 13.81 □ 

Theorem 4.5. The function Mb(A/") defined above on ported unimodular oriented ma- 
troids satisfies the properties: (Symbols like E and P denote sequences and M^A/") de- 
pends on the g e , r e and e.) 

1. If e G E is neither a separator nor a port, and E' = E\ E, then 

M E (Af) = e(PE)e(PE') (g e M E > (A/ - / e) + r e M E >(J\f \ e)) . 
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2. If Mi{P\i Ei) and N' 2 (P 2 , E 2 ) have disjoint ground sets and E — E\ UE 2 , then 

M £ (M®M) = e(P 1 P 2 E)e(P 1 E 1 )e(P 2 E 2 ) M El (M) Mj^M). 
5. If P = $ and B(N) denotes the collection of bases of N, then 

M E {M) = e{E) 9bTb. 

BeB(Af) 

Proof. Immediate from Theorem 13.31 and the above remarks. □ 

5 Basis, Set and Flat Expansions 

Theorem 14 . 51 shows that when Af(P, E) is a unimodular matroid, M^JV) is a substitution 
of extensors and u — v — into Rp(J\f) of Definition 14.11 Proposition 14.21 demonstrates 
Rp is a ported Tutte function. This motivates general study of ported and parametrized 
Tutte functions of matroids. Below, we extend some known general expressions for Tutte 
functions and invariants. Our results are more clearly expressed and no harder to prove 
than if the parameters were omitted. Furthermore, we identify an unsolved problem due 
to combining parameters with ports. 

Evidently, each Tutte function value F{M) is determined via the Tutte equations from 
parameter values and from values of F on decomposable minors of N '. One way to study 
a Tutte function is to present the definition independently of the Tutte equations and 
then prove that the function so defined satisfies the Tutte equations. A second way is to 
specify the parameters and the function values on indecomposable minors and then prove 
a solution exists and is unique. In other words, one proves that all the Tutte equations 
are consistent with the given parameters and values, and that the solution is unique. 
(This extends to our Tutte functions remarks of Pak [34].) The issue that arbitrary 
parameters and values on indecomposable sometimes fail to be consistent or sometimes 
result in multiple solutions was studied in [5,23,67]. Fortunately, the Tutte functions (of 
matroids, oriented matroids and of extensors) studied in this paper are all defined the 
first way. However, our work suggests the open problem to extend these studies to the 
ported Tutte equations. The indecomposables will then include connected matroids or 
oriented matroids with more than one element, i.e., not just loops and coloops. 

Zaslavsky [67] defines the normal class of parametrized (but not ported) Tutte func- 
tions as those Tutte functions for which there exist u and v for which, for all e, the point 
value on coloop e is r e u + g e and the point value on loop e is g e v + r e . The normal Tutte 
functions are exactly those obtained by substitutions into the parametrized corank-nullity 
polynomial. All Tutte invariants are normal Tutte functions and non-normal Tutte func- 
tions do not express much of the matroid structure — See [5,23,67] for details about how 
parametrization complicates Tutte invariant theory. 

Our Theorem 14.41 expresses how M^TV) fits into the natural ported generalization of 
the normal class. Our extensor Tutte function of ported oriented unimodular matroids and 
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its invariant specialization is expressed by a substitution into the ported corank-nullity 
polynomial of oriented matroids. 

While only the normal Tutte functions have corank-nullity polynomial expressions, 
they all have basis expansion expressions [67]. In the rest of this section, we discuss these 
and other expansion expressions for ported unoriented and oriented matroids. 

The basis expansion originated by Tutte [51] for graphs and Crapo [21] for matroids 
depends on a particular but arbitrary ground set element order O. Each basis determines 
a term from the internal and external activities of elements with respect to that basis 
according to the ordering O. Our way to generalize is to restrict O to orders in which 
every port element is ordered before each e G E. (We use the convention that the 
deleted/contracted element is the last, i.e., greatest element under order O eligible for 
reduction.) 

Gordan and McMahon define [26] a "computation tree" to formalize the application of 
a subset of Tutte equations to a matroid and some of its minors. Each (Tutte) computation 
tree determines a polynomial in the parameters and point values. Therefore, when M is 
in the domain of a Tutte function, each of these computation trees determine the same 
value. Computation trees are a way to give a basis expansion expression in terms of a more 
general definition of internal and external activities of elements with respect to a basis. 
The expansion is more general because it is based on any Tutte equation computation 
rather than on an element order O. We will extend to ported computation trees the 
classification [26] of elements as internally or externally, active or passive with respect to 
each path down the tree. In each case, the result is an interval partition of the boolean 
subset lattice of E. 

Definition 5.1. Given J\f(P,E), a P-subbasis F G Bp(N) is an independent set with 
F C E (so F PI P = 0) for which F U P is a spanning set for M{P, E) (in other words, F 
spans N '/P, see [37].) 

Proposition 5.2. For every P-subbasis F there exists an independent set Q C P that 
extends F to a basis F UQ G B(M). Conversely, if B e B(M) then F = B n E = B \ P 
is a P-subbasis. 

Proof. Immediate. □ 

Definition 5.3 (Activities with respect to a P-subbasis and an element ordering 

O). Let ordering O have every p G P before every e G E. Let F be a P-subbasis. Let B 
be any basis for M with F C B. 

• Element e G F is internally active if e is the least element within its principal cocir- 
cuit with respect to B. Thus, this principal cocircuit contains no ports. The reader 
can verify this definition is independent of the B chosen to extend F. Elements 
e G F that are not internally active are called internally inactive. 

• Dually, element e G E with e ^ F is externally active if e is the least element 
within its principal circuit with respect to B. Thus, each externally active element 
is spanned by F. Elements e G E \ F that are not externally active are called 
externally inactive. 
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Definition 5.4 (Computation Tree, following [26]). A ported (Tutte) computation 
tree for Af(P, P) is a binary tree whose root is labeled by Af and which satisfies: 

1. If Af has non-separating elements not in P, then the root has two subtrees and there 
exists one such element e for which one subtree is a computation tree for Af /e and 
the other subtree is a computation tree for Af\e. 

The branch to Af je is labeled with "e contracted" and the other branch is labeled 
"e deleted". 

2. Otherwise (i.e., every element in S(Af) \ P) is separating) the root is a leaf. 

An immediate consequence is 

Proposition 5.5. Each leaf of a P-ported computation tree for Af(P, E) is labeled by the 
direct sum of some minor of Af on P (oriented if Af is oriented) summed with loop and/or 
coloop matroids with ground sets {e} for various distinct e G E (possibly none). 

Definition 5.6 (Activities with respect to a leaf). For a ported computation tree 
for Af(P, E), a given leaf, and the path from the root to this leaf: 

• Each e G E labeled "contracted" along this path is called internally passive. 

• Each coloop e G E in the leaf's matroid is called internally active. 

• Each e G E labeled "deleted" along this path is called externally passive. 

• Each loop e G E in the leaf's matroid is called externally active. 

Proposition 5.7. Given a leaf of a ported computation tree for Af(P,E): The set of 
internally active or internally passive elements constitute a P-subbasis of Af which we say 
belongs to the leaf. Furthermore, every P-subbasis F of M belongs to a unique leaf. 

Proof. For the purpose of this proof, let us extend Definition 15.61 so that, given a compu- 
tation tree with a given node % labeled by matroid Afi, e G E is called internally passive 
when e is labeled "contracted" along the path from root Af to node i. Let IPi denote the 
set of such internally passive elements. 

It is easy to prove by induction on the length of the root to node i path that (1) 
IPi U S(J\fi) spans Af and (2) IPi is an independent set in Af. The proof of (1) uses the 
fact that elements labeled deleted are non-separators. The proof of (2) uses the fact that 
for each non-separator / G Af/IPi, f U IPi is independent in Af. 

These properties applied to a leaf demonstrate the first conclusion, since each e G E 
in the leaf's matroid must be a separator by Definition 15.41 

Given a P-subbasis P, we can find the unique leaf as follows: Beginning at the root, 
descend the tree according to the rule: At each branch node, descend along the edge 
labeled "e-contracted" if e G P and along the edge labeled "e-deleted" otherwise (when 
e^F). (This algorithm also operates on arbitrary F' C P.) □ 
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The above definitions and properties enable us to conclude: 

Proposition 5.8. Given element ordering O in which every p G P is ordered before each 
e $l P , suppose we construct the unique P-ported computation tree T in which the greatest 
non-separator e G E is deleted and contracted in the matroid of each tree node. 

The activity of each e G E relative to ordering O and P-subbasis F C E is the same 
as the activity of e defined with respect to the leaf belonging to F in T. 

Definition 5.9. Given a computation tree for (oriented) matroid Af(P,E), each P- 
subbasis F C E is associated with the following subsets of non-port elements defined 
according to Definition 15 . 61 from the unique leaf determined by the algorithm given above. 

• IA(F) C F denotes the set of internally active elements, 

• IP(F) C F denotes the set of internally passive elements, 

• EA(F) C E\F denotes the set of externally active elements, and 

• EP(F) C E \ F denotes the set of externally passive elements. 

Proposition 5.10. Given a computation tree for W(P, E) , the boolean lattice of subsets 
of E is partitioned by the collection of intervals [IP(F), F U EA(F)] (note FUEA(F) = 
IP(F) U A(F)) determined from the collection of P-subbases F, which correspond to the 
leaves. 

Proof. Every subset F' C E = S(Af) \ P belongs to the unique interval corresponding to 
the unique leaf found by the tree descending algorithm given at the end of the previous 
proof. □ 

Dualizing, we obtain: 

Proposition 5.11. Given a computation tree for J\f(P, E) , the boolean lattice of subsets 
of E is also partitioned by the collection of intervals [EP(F), E \ F U IA(F)] (note E \ 
F U IA(F) = EP(F) U A(F) ). 

Proof. The dual of the tree descending algorithm is to descend along the edge labeled 



The following generalizes the basis expansion expression given in [67] to ported (ori- 
ented) matroids, as well as Theorem 8.1 of [37]. 

Definition 5.12. Given parameters g e , r e , point values x e , y e , and (oriented) J\f(P, E) the 
Tutte polynomial expression determined by the sets in Definition 15 . 91 from a computation 
tree is equal to 



Each Tutte polynomial expression is constructed by applying some of the Tutte equa- 
tions. Therefore, if Af(P, E) is in the domain of Tutte function /, then f{M) is given by 
any Tutte polynomial expression with f(N/F\P) substituted for each oriented or unori- 
ented matroid monomial [M / F\P]. (This generalizes the expression used in [67] to define 
the Tutte polynomial under the condition that all expansions yield the same expression.) 



e-deleted" if e G F. 



□ 




(23) 



F&Bp 
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5.1 Boolean Interval Expansion 



The following proposition expresses the ported corank-nullity polynomial in terms of a P- 
subbasis expansion. It is obtained by substituting binomials x e = g e +r e u, y e = r e +g e v and 
leaving the matroid variables unchanged in Definition l5.12l The different expansions from 
different element orderings and Tutte computation trees all express the same polynomial 
because Proposition 14.21 demonstrates that Rp is a ported Tutte function and the values 
of Rp on coloop, loop and indecomposable matroids are readily verified to be given by 
these substitutions. 

Proposition 5.13. The polynomial Rp(N) is given by the following activities and boolean 
interval expansion formula: 

RpW = E W/F\P\ ( E 9ku(e\f\l) r Lu{F \ K) ) (24) 

F&Bp IP(F)CKCF 
EP(F)CLCE\F 

Proof. Let A — KU(E\F\L) within the above expansion. We can verify A — E \ A — 
LU (F\K). For each A C E a unique P-subbasis F, and two tree leaves are determined, 
one by the tree descending algorithm and the other leaf by the dual algorithm. Thus A 
and A respectively belong to intervals within the boolean lattice partitions of Propositions 
EM and ETTTI In particular, A e [IP(F) , F U EA(F)} and A e [EP(F) ,E\FUIA(F)]. 
Therefore the terms in the above sum are equal one by one to the terms in the corank- 
nullity polynomial's subset expansion (Definition 14. □ 

For the purposes of this paper it was sufficient to recognize that our extensor valued 
ported parametrized Tutte function of unimodular oriented matroids belongs to the nat- 
ural generalization of Zaslavsky's normal class. As such, it has, for arbitrary parameters, 
expressions obtained by substitutions into (1) computation trees, (2) ported parametrized 
Tutte polynomials from such trees, (3) various P-subbasis expansions, and (4) the ported 
parametrized corank-nullity polynomial. 

We state here the open problem to include ports into the results of Zaslavsky, Bollobas 
and Riordan, and Ellis-Monaghan and Traldi: Can we classify with universal forms all of 
the ported parametrized Tutte functions according to their parameters, non-port point 
values, and the values on oriented or unoriented minors on port ground sets? 

5.2 Geometric Lattice Flat Expansion 

A formula for the unparametrized ported Tutte (or corank-nullity) polynomials of non- 
oriented matroids in terms of the lattice of flats (closed sets) and its Mobius function 
was given in [12]. We generalize: (1) The expansion's monomials [Q] can signify either 
oriented matroid minors, when M is oriented, or non-oriented minors when M is not 
oriented. (2) The polynomial is parametrized with r e , g e for each e G E. The derivation 
relies on the fact that the oriented or non-oriented matroid minor [J\f/A\P] (according to 
whether M is oriented or not) depends only on the flat spanned by A C E. 



33 



Proposition 5.14. Let J\f(P, E) be an oriented or unoriented. Let Rp(N) be given from 
Definition \4-l\ In the formula below, F and G range over the geometric lattice of flats 
contained in E. 

R P (Af)(u, v) = Y}Q] E u^-o Q -» F v-i> F E MG, F) ]J(r e + g e v) (25) 

Q F<E G<F eSG 

W/F\P] = [Q] 

Proof. It follows the steps for theorem 8 in [12]. □ 

Remark: The chirotope values for the oriented matroid minor N' = M / F\P are 
Xj\f'(X) = xn(XB f ) where X is restricted to sequences over P and B F is any basis 
for the fiat spanned by F. While this formula defines a chirotope function only up 
to a constant sign factor, the oriented matroid (which is what the monomial [M / F\P] 
denotes) is uniquely defined. We mention this because when we evaluate the corank- 
nullity polynomial to obtain M^AT) for the unimodular oriented matroid M substitute 
an extensor for each [Af/A\P\. However, the object we substitute is M0(N/A|P), not 
N/A\P. It is the unique extensor defined by equation (jlHj) applied to one of the chirotopes 
that present [7V/A|P] (or to any other representation of [A/"/A|P] for that matter). We 
already remarked that equation (fT^j) is unchanged when its argument changes sign! 

6 Electrical Networks 

The most common way to formulate the electrical network problem that occurs in our 
subject is to use the discrete and parametrized Laplace's equation. The insights from 
matroid theory appear with greater clarity when the analyses of the solution, which is 
in terms of the electrical potential values at graph vertices, are replaced by analyses in 
terms of voltage and current variables directly associated with edges. After we review the 
Laplacian, we move on to edge-based formulations. Some details of relationships between 
our results and electrical network problems are then demonstrated. 

The following observations underlie the combination of ideas about electrical networks 
and ported oriented matroids: First, the basis enumerator is a Tutte function that hap- 
pens to be a determinant in the case of unimodular (or regular) matroids, such as the 
graphic matroids. Second, topics involving determinants, including the chirotope presen- 
tation of realizable oriented matroids, can profitably be studied with exterior algebra. 
Third, non-trivial electrical networks (see ^6.2|) must have port elements for supplying 
power or for connecting to an external environment, in addition to the resistor elements 
usually modeled by graph edges. We are interested in their combinatorial properties be- 
yond spanning tree counts [9-15]. Finally, electrical flows and potential differences are 
inherently directional. The patterns of their directions (expressed by sign functions on 
graph edges) feasible under Kirchhoff's two laws are precisely the vector and covector 
families, respectively, of the graphic oriented matroid. Indeed, the duality between these 
laws, current and voltage, is characterized by oriented matroid theory. 
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6.1 Discrete Laplacian 



The combinatorial (or discrete) Laplacian is the matrix of coefficients in the equations 
(|2T)j) below in variables 1 < i < n. These discrete Laplace equations model (among 
other situations) a resistive electrical network when <f>i represents the electrical potential 
(or voltage) at vertex i and constant Ii represents the current flowing into vertex i. 



If each conductance g e is non-negative, or is either zero or generic, then the rank of 
the n x n Laplacian matrix is n — k, where k is the number of connected components in 
the n vertex undirected graph whose edges are the e = ij with g e ^ 0. Each order n — k 
non-singular diagonal submatrix is called a reduced Laplacian. (In ^6.21 we will express 
Laplace's equations with an M^(N) independent of k.) The reduced Laplace equations, 
together with 0i = for each vertex i corresponding to a deleted column, model a network 
where each such i is contracted into a single grounded vertex whose potential is fixed to 
zero and whose external current is unrestricted 3 . Since the graph is undirected, equations 
(I26|) imply that the current into the grounded vertex equals the sum of the Ij for the 
non-grounded vertices. 

The inverse of a reduced Laplacian matrix is called the discrete Green's function in [20]. 
This inverse matrix's elements are each expressed (using Cramer's rule) by a ratio of an 
order n — k — 1 minor to a common order n — k minor denominator. The list of all minors, 
of all orders, is an example of Pliicker coordinates — Here, these are the ( n ) maximal 
minors of the matrix obtained by appending the n x n identity matrix to the side of the 
Laplacian. 

The Matrix Tree Theorem asserts that each n — 1 order minor equals ±Y2dT, the 
enumeration of spanning trees T by products of edge parameters gr = FleeT^e- See [9] 
for similar interpretations of all the minors and for generalizations to directed graphs. The 
formulas we call "Maxwell's rules" were given without proof by Maxwell [41], and the dual 
forms of them were proved by Kirchhoff [33]. Maxwell also described the static equilibrium 
solution for stressed linear elastic framework in terms of enumerations over minimally rigid 
subframeworks [42]; this enumerated set is the basis set for the rigidity matroid [27]. The 
one-dimensional case is analogous to the electrical problem. The survey by Biggs [3] 
covers the discrete Laplacian, the Matrix Tree Theorem, and the use of spanning tree 
enumeration to solve the discrete Laplace equations, and many additional topics, including 
the asymmetric discrete Laplacian. Biggs presents the Kirchhoff's solution method, and 
Nerode and Shank [44], also used by Bott and Duffin [6], Smith [47] and Maurer [40]. 
This method constructs a symmetric projection matrix from a sum of fundamental cocycle 
matrices, one for each spanning tree. Analysis of basis exchange, i.e., the pivot calculation 
implies the appropriately weighted matrix sum is symmetric. We plan to present the 
generalization of this argument to extensors in a future publication. 

3 Do not confuse with ground set. 




(26) 



{j\e=ij£E} 
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Tree counting, the discrete Laplacian and electrical network models with parameters 
have a spectrum of applications including electrical circuit theory, knot theory, random 
walks and the analysis of Markov chains (see for example [3,19,20,22]). Their application 
to square dissections is described in [7, 53] ; Tutte gives a Laplacian based "barycentric 
embedding" proof of Kuratowski's Theorem in [52]. 

It is generally known among electrical engineers in circuit theory that the same kinds of 
homogeneous rational polynomial functions that appear in Maxwell's rule occur generally 
as the coefficients (and minors of them) in all of the linear relationships between the port 
quantities that define the externally observable characteristics of a linear resistive network. 
Our results display this principle within the mathematical contexts of the enumerative 
combinatorics of graphs, oriented matroids and exterior algebra: See Corollary 13. 71 Some 
electrical network analysis software actually enumerates trees and related structures to 
do "symbolic analysis." See for example [14,24,49,50]. 

6.2 Electrical Network Equations in terms of Edges 

Consider a graph with two kinds of edges, called ports P and resistors E. The graph is 
directed with an arbitrary edge orientation. Let unimodular extensor N(P, E) present its 
ported graphic oriented matroid. Let r be the rank of this matroid. 

Let g e ,r e be parameters for each e G E. The extensors i(N)(P t UP„ UP) and 
w(N J -)(P t UP„ UP) defined in ^determine the electrical network equations in the way 
expressed by Theorem 12.31 applied to N -1 . 

To be specific, these equations are a linear system on the \E\ + 2|P| variables {x e , . . . ; 
v p , . . . ; i p , . . .}. Each e G E is associated to variable x e . Each p G P is associated to two 
variables, v p called the voltage and i p called the current. Let matrix K with r rows be 
any matrix that presents t(N). Let C be any matrix with |P| + |P| — r rows that presents 
t^N -1 ). These matrices express Kirchhoff's equations combined with a homogeneous 
expression of Ohm's law. K determines the following current equations: 

^2 K 3^P + S K 3fi X * = for J = 1, . . . , T. 

pGP eG-E 

C determines the following voltage equations: 

^ Cj,p v v p + ^2 Cj, e x e = for j = 1, . . . , \E\ + |P| - r. 

peP e£E 

It is helpful to see the electrical network equations in terms of N directly. Let N be any 
matrix presentation of N; for example, a reduced oriented incidence matrix of the graph. 
Let N x present N -1 ; the rows of A r± comprise a basis for the cycle space of the graph. 
The current equations can be written: 

^2 N 3,P i P + ^2 N h^eX e = for j = 1, . . . , V. 

pgP eG-B 
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The voltage equations can be written: 

Yl N t v v + Yl N t* r * x * = for j = 1, . . . , \E\ + \P\ - r. 

peP e&E 

The equations which M^(N) presents are obtained by eliminating all the variables x e , 
e G E in the voltage and current equations taken together. Corollary 13.61 indicates the 
rank of the resulting system of \P\ equations on 2\P\ variables is \P\, provided that the 
parameters are generic or all positive. 

The above analysis illustrates the role for the port element distinction in modeling a 
physical system. Each non-port element models a completely defined subsystem. The 
"proto- voltage [47]" x e parametrizes the state of one electrical resistor, for example. The 
behavior of this resistor is thus defined by Ohm's law: When g e and r e are both non-zero, 
the current is g e x e if and only if the voltage is r e x e . The entire model (the graph, for 
example) specifies all the interactions (via Kirchhoff's laws, for electricity) between its 
subsystems. Each port element models an interface pertaining to an interaction of the 
system with an unspecified environment, for observing the system behavior of interest 
to the application, and to help specify how certain larger systems are composed out 
of previously entire subsystems. For us, the environment is assumed, for each port, 
to constrain the currents into one terminal and out of the other terminal to be equal. 
Environmental constraints between voltages at terminals belonging to distinct ports are 
forbidden as well. (Engineering models encompass multiport elements, whose behavior 
is specified using multiple port elements [16,18,46]. For example, a linear multiport 
element is specified a linear constraints among the variables associated with its ports; this 
generalizes Ohm's law to so-called multi-terminal resistors. Each of our ported objects 
can model a single multiport element within a larger model. A topic for future research 
is to abstract this along the lines given here.) 

Let a graph on vertices {1, . . . ,n} be given with conductance parameters g e for each 
edge. We now derive Laplace's equations from the voltage and current equations. Let us 
append a new vertex (which will be grounded) and n port edges pi G P, with each p^ 
directed from vertex to vertex i. For simplicity take parameter r e = 1 for each edge e 
from the original graph. We can choose 7V so that the current equations are 

^ ^ Jp,e9e-Eei P £ P 

and the voltage equations are 

n 

x e = <P Pr - 4> P s = Y ^ViA-Pii where e = rs, e G E. 

i=l 

where we used potential (relative to vertex 0) 0j in place of v Pi and J is the oriented vertex- 
edge incidence matrix of the original graph. Laplace's equation is obtained by eliminating 
the variables x e which represent differences of potential across resistor edges. Indeed, one 
presentation for M^N) in this case is the n x 2n matrix [I re A] formed by concatenating 
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the identity matrix with the Laplacian matrix A. One manifestation of Theorem 13.31 is 
therefore that each of the forest enumerating polynomials given by an arbitrary minor, of 
any order, of the Laplacian is a (non-ported) Tutte function of graphic oriented matroids. 

6.3 Maxwell's Rules 

Kirchhoff [33] first described the solution of a resistive electrical network problem in 
terms of enumerations of spanning trees and certain forests (or dual forests). One such 
description is the following surprising yet classical result, one of the "easily remembered" 
rules stated without proof by Maxwell [41]. 

Let's discuss Maxwell's rule for a network with one port edge p in terms of graph edges. 
Let all edges ^ p represent unit resistors in an electrical network. Port edge p is to demarks 
two terminal vertices. Maxwell's rule asserts that the equivalent resistance between the 
two terminal vertices equals the quotient of the count A of spanning trees that contain 
edge p divided by the count B of spanning trees that omit edge p. Note that whenever 
a two-tree spanning forest F for which F Up is a tree counted by A is contracted and 
the other non-p edges deleted, p becomes a coloop. When this is done for B, p becomes 
a loop. A and B are the two Pliicker coordinates of the extensor we construct in this 
paper. They are the coefficients in the linear equation Ai p + B(—v p ) = that relates 
the port current and the (negated) port voltage. The reader can verify that each of A 
and B satisfy the additive ported Tutte equation; so (A, B) satisfies it. (It is instructive 
to calculate how (A,B) behaves under the multiplicative ported Tutte equation.) The 
spanning tree count resistance expression (A, B) generalizes with weighted spanning tree 
generating function when conductance parameters g e are given for e ^ p. 

Suppose two ports, p\ and P2 are given. One, say p±, denotes the terminal pair 
between which a specified current value i Pl is constrained to flow (assuming this is feasible.) 
The other, p 2 , denotes the terminals between which to observe or measure the potential 
difference v p2 (voltage) that results when the specified current constraint is the only 
electrical power source (assuming this voltage value is unique.) Maxwell's rule for one 
port is the specialization with pi = P2 = P- 

The two-port Maxwell's rule holds that the port variables are related by Ai Pl + 
B(—v P2 ) = were A and B are again generating functions, except that some of the 
monomials in A might be negative. The spanning trees that omit p\ and p 2 are enumer- 
ated by B. Term A enumerates forests F for which each of FUpi and FUp 2 is a spanning 
tree are enumerated by A. For each such F, FU {pi,P2} contains a unique circuit, which 
of course contains both p\ and p 2 . The relative directions in which p\ and p 2 are traversed 
determines the sign of rkgp in A. 

It is amusing go back to the one port p case and derive the Maxwell's rule coefficients 
(A, B) from the principles we just illustrated. The port behavior of the network consisting 
of one coloop p is defined by 1 • % p + • v p = 0: This constraint can be expressed by Pliicker 
coordinates (1,0). The port current is and the port voltage is unconstrained since the 
graph has no cycles. Dually, the port behavior of the one port loop network is expressed 
by (0,1). The port voltage is but the current is unconstrained. One consequence of 
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our main theorem proves that A(l, 0) + 5(0,1) = (A, B) gives the Pliicker coordinates 
of the constraint that the original graph imposes on its port current and (negated) port 
voltage variables, in accordance with Maxwell's rule Ai p + B(—v p ) = 0. Our theory 
justifies the choices of (1,0) and (0,1) rather than any other multiple (a, 0) or (0, 0) in 
A(a,0) + B(0,p). 

6.4 Deriving Maxwell's Rule 

See [16] for an elementary derivation that includes the version of Maxwell's rule [41] that 
applies to to a pair of ports that do not share a common vertex. Graphic matroid orien- 
tation becomes relevant in this situation. Explicit port edges have proven their usefulness 
in electrical network analysis [18, §13.6]. In this situation, some Pliicker coordinates of 
M.e(N), as polynomials in g e , r e , have terms of opposite sign only when two port edges do 
not share a vertex. Our contribution is to characterize the signs within the theory of ori- 
ented matroids and Tutte functions. It is true that such polynomials can be expressed in 
terms of minors of the Laplacian; this was done by manipulation of solutions to Laplace's 
equation in [7]; see also [53]. However, our extensor and oriented matroid formulation 
enables the analysis to be done without the introduction of vertices. 

Our derivation of Maxwell's rule for two ports (of which the one port version is a special 
case) begins with the electrical network equations with P = {pi,P2}- Let M = M^(N) 
as defined in ^be as discussed above, and let M be any 2x4 matrix presentation of M. 
The two equations 

k 

M 12 =0 

. V ? . 

are obtained by eliminating the variables x e , e G E from the electrical network equations. 
The currents i\, %2 in edges p\ = ah and P2 = cd flow from vertices a to b, and c to d 
respectively. The voltage (drop) v\ across edge pi is the potential at a minus the potential 
at b; the corresponding convention defines the voltage v% across edge p2- 

We will assume that all the r e — 1 and that M[p ll) p 2 „] ^ 0. The latter is assured 
from Corollary 13.101 (4.) provided that E contains a spanning tree and all the g e are 
either positive or generic. Under these conditions, the transfer resistance pi\ given by 
{—V2)l%\ when %i = and i\ ^ is well-defined. (These conventions are used so that when 
Pi and P2 are identical or parallel, p signifies the familiar equivalent resistance which is 
always positive when E is connected and all g e > 0.) 

Proposition 6.1 (Maxwell's Rule). Given the electrical network graph model described 
above, let B denote the collection of edge sets T C E of trees that span the vertex set V , 
and assume J2 T&B g T ^ 0. 

For vertices i,j,k,l, let Bikji be the collection of all F C E for which the subgraph 
(V, F) is a forest with exactly two trees where vertices i and k are in one tree and j and 
I are in the other tree. 
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The transfer resistance p 2 \, where p\ = ab and p 2 = cd, is well-defined and is given 

by: 



T^FeB acM 9F Y.FeB ad . bc 9f 

P21 = : • (27) 

Proof. We will abuse the notation slightly by using Vk and ik for ground set elements pk v 
and pkt,, and Vfc, i& for the corresponding extensors, k = 1,2. 

Corollary 13 .1 01 ( 4.) shows that M[t>it> 2 ] (= Nl[pi v p 2v \) = Y1t&b9 t ^ ^' so by Cramer's 
rule, 

_V2 ( M g (N)[mii] 
P21 ~ ii V M E (N)[ Vl v 2 } 
Let us apply Corollary 13.81 to the numerator and denominator. To do this, we first 
calculate M0(N) for 6 unimodular extensors N({pi,p 2 },$) that present the 6 oriented 
matroids on P = {pi,p 2 }, which are all graphic. For each of the 6, we can then determine 
the Pliicker coordinate values M0(N)[t> 1 i 1 ] and M0(N)[t>iU 2 ]. Those oriented matroid 
minors for which one of these values is non-zero will characterize, together with the 
rank conditions in Corollary 13.81 which forests or trees contribute to each sum. These 
characterizations of the forest or tree terms, and of their signs, analyzed for each Pliicker 
coordinate, will complete the derivation of Maxwell's rule. 

Four of the 6 oriented matroids are the direct sums of either the loop A/oQoi) or coloop 
Miijpi) with either the loop or coloop on p 2 . The other two oriented matroids are the ori- 
entations of the 2-circuit matroid on {pi,P2}- Let Mi denote the oriented circuit ±(H — ); 
graphically, p\ and p 2 are parallel. So denotes the oriented matroid of antiparallel 
Pi and p 2 . Table [T] lists these six distinct ported oriented matroids J\f{{pi,p 2 }, 0), their 
unimodular extensor presentations N({pi,p 2 }), and the corresponding extensor values 
M0(N). The M0(N) values are easily found up to sign. The signs are given in Proposi- 



Table 1: The six oriented matroids on {pi,p 2 }- 



matroid 


N 


M (N) 


A/" (pi)eA/'o(p 2 ) 


±1 




viv 2 


A/" (pi) ©M(P?) 


±P2 




I2V1 


M(Pi) ©A/" (p 2 ) 


±Pl 




iiv 2 


M(p0 ©M(p 2 ) 


±PlP2 




iii2 


A/f 


±(Pl - P2) 


(ii - i 2 )(vi + v 2 ) 


= iiVi + i x v 2 - i 2 vi - i 2 v 2 




±(Pl +P 2 ) 


(ii + i 2 )(v 2 - Vi) 


= iiv 2 - iiVi + i 2 v 2 - i 2 vi 



tion 13.51 One way to calculate is to analyze the corresponding electrical network with 2 
ports and no resistors. For example, the network with oriented matroid Ao(pi) ® N\(p 2 ) 
constrains its port p±, a loop, to have voltage drop v 1 = but its port p 2 , a coloop, to 
have current i 2 = 0. The current in the loop and voltage across the coloop are uncon- 
strained. The solution subspace corresponds to equations (t>i = 0; i 2 = 0) . The extensors 
representing these equations are avii 2 , a^O. 
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Similarly, the network of two parallel ports (case J\f^~) constrains the sum of voltage 
drops going around the oriented circuit to be 0, so Kirchhoff's voltage law is expressed by 
v \ — v 2 = 0. Kirchhoff's current law in the same network is expressed z'i + i 2 = 0. Hence 
the corresponding extensor is ±(vi — v 2 )(ii + 12)- 

We complete the derivation. First for the denominator. From tabled the only terms 
in (}2"T]) of Corollary 13. 81 that might contribute to M^N)^]/^] are those for which M / A\P 
is the matroid of two loops Ao(pi) © Nq{P2) because the only appearance of vxv 2 is in 
that matroid's row. The rank conditions further restrict the contributing A to spanning 
trees. 

Finally, for the numerator M^(N)[uiii], we locate ±Viii in the bottom two rows. 
These appearances have opposite sign. For ACE with M /A\P = J\f[~, the contribution 
is —gA- The sign is opposite when J\f/A\P = A/"] 4 ", so the distinct orientations of the 
2-circuit obtained when contracting A account for the opposite signs in (j27jl . We can 
again verify from the rank conditions that the F = A contributing to the numerator of 
(|27jl are the spanning forests with 2 trees containing the indicted vertices as claimed. 

Note that the sign dependence of Mb(N) on e and the order of P = p\p 2 cancels in 
the ratio p 2 i- □ 

While Theorem 16 .11 can be proved by elementary arguments as in [16], the above proof 
demonstrates how it can be derived from the forgoing theory using algebraic calculations. 

Remark: The one port version is immediately derived using a graph where p\ and p 2 
are parallel edges because our proof puts no special conditions on the two ports. 

6.5 Signed Contributions of Spanning Forests 

In some cases of graphs given with two or more port edges, some of these Pliicker coordi- 
nates will equal the difference between the counts of two kinds of spanning forests. Such 
a coordinate pertains to the coefficient that relates a quantity observed at one port to of 
a voltage or a current quantity at a different port. Our results show how oriented ma- 
troid properties determine the sign by which each forest F contributes to this coordinate 
(Corollary 13. 8j) . In particular, that sign is determined by the graphic oriented matroid on 
the (directed) port edges obtained by contracting F and deleting the remaining non-port 
edges. The simplest case where distinct signs do occur is when both orientations of the 
same 2-circuit matroid on two ports appear in this process. The contribution of F to the 
Pliicker coordinate with a given index X is calculated in a particularly simple way: We 
solve the electrical network of port edges only (with no resistors!) that resulted from this 
deletion and contraction, after checking matroid rank conditions necessary for at least one 
coordinate to be non-zero. The contribution equals the Pliicker coordinate with the same 
index X from the solution of the latter, port-edge-only network, weighted by gF^E\F- It 
should be noted that the sign of the contribution is determined by the oriented graphic 
matroid on ports only, independent of the particular F contracted to obtain this oriented 
matroid. Details appeared in our proof of Maxwell's rule for two ports, section 16.31 
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7 Additional Background and Directions 



7.1 Ground Set Orientation 

The ground set orientation and its role in defining a canonical dual of an extensor, and 
our M^(N) are motivated by the idea of orientations of orientable manifolds and the 
definition of pseudo- forms (or "forms of odd- kind" attributed to de Rham in [25]) in 
the mathematics of physics. A pseudo-form is an antisymmetric multilinear operator 
f = f e that is parametrized by the orientation e and for which ef e is independent of 
the orientation [25]. So, ef e is a well-defined form. In physics, an orientation specifies 
one's convention, say by a right-handed coordinate system, for how one defines a positive 
volume or other naturally unsigned physical quantity in terms of an exterior algebra form. 
The orientation specifies which ordered bases determine right handed coordinate systems. 
In this context, the orientation e is a ±1 function for which e(-Bi)e(-B 2 ) is the sign of the 
determinant of the local Jacobian matrix which relates the ordered bases Bi,B 2 . 

7.2 Computational Complexity 

Among the non-trivial Tutte invariant functions of succinctly presented graphs or ma- 
troids, only two (unless V = #V) are polynomial time computable [32,57]. One such 
function, the number of bases, is computable by the Matrix Tree Theorem for graphs 
and its extension to unimodular matroids. This number is well-known as the evalua- 
tion T(Af, 1,1) of the Tutte polynomial function T(Af,x,y) of matroids M. However, 
computing T(J\f, 1, 1) is #P-complete for arbitrary non-unimodular matroids [56]; this 
follows because counting the perfect matchings in a bipartite graph is a #P-complete 
problem [55]. 

The other easy-to-compute invariant is determined by the dimension of the intersection 
of a linear subspace and its orthogonal complement over a finite field [57]. The cited 
papers prove that all of the other Tutte matroid invariants are either trivial or #V- 
complete. More recently, analogous computational hardness results have been proven 
for Tutte functions of graphs (thus implying their hardness for matroids). Among these 
results, is that evaluating the parametrized Tutte polynomial for given matroids is a VNP- 
complete problem [38]. Here, Valiant's non-uniform algebraic complexity model [54] is 
used, which counts as one deterministic step each evaluation of a polynomial on constants, 
variables or previously computed values. VNP is this model's class that is analogous to 
NP in the Turing machine model. (See the references in [38]). 

A full account of the computational complexity of Tutte invariants of graphs and 
matroids is given in [32,45,57,59-61]. 

We remark for computationally-inclined readers that: 

1. The Tutte equations describe non- unique recursive algorithms to compute Tutte 
functions that generally require 2'^' steps. 

2. A \P\ x 2\P\ matrix representing our extensor can be computed from a graph or 
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locally or totally unimodular matrix presentation of a ported oriented unimodular 
matroid. 

One suitable algorithm is simple matrix block manipulations followed by Gaussian 
elimination. Such elimination-based algorithms use polynomial bounded numbers 
of field operations. Therefore, computation of our extensor generalization of the 
basis enumerator on graphic and other unimodular matroids is a polynomial time 
problem when all r e = g e = 1. 

7.3 Other Directions 

It is natural to generalize the current and voltage equations so their respective solutions 
subspaces (taken to be within KS) are not orthogonal [10]. This leads to the directed 
graph version of the Matrix Tree Theorem. It did lead as well to a "oriented matroid pair" 
model for combinatorial conditions for certain equations with monotone non-linearities to 
be uniquely solvable [13]. These conditions were stated in terms of two oriented ma- 
troids with a common ground set having complementary rank and no common non-zero 
covector; the current paper provides the insight that these two were obtained by dele- 
tion/contractions to eliminate port elements. Investigations of a generalization of the 
Tutte polynomial to two matroids with a common ground set were also begun in [62]. 

The computation tree formalism was used in greedoid generalizations [26]. of the Tutte 
polynomial because those generalizations do not always have an activities expansions 
based on element orders. We leave investigation of "ported greedoids" to the future. 
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